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Abstract 

Numerous aspects and mechanisms of color confinement in QCD 
are surveyed. After a gauge-invariant definition of order parameters, 
the phenomenon is formulated in the language of field correlators, to 
select a particular correlator responsible for confinement. In terms of 
effective Lagrangians confinement is viewed upon as a dual Meissner 
reffect and a quantitative correspondence is established via the popu- 
lar abelian projection method, which is explained in detail. To deter- 
mine the field configurations possibly responsible for confinement, the 
search is made among the classical solutions, and selection criterium 
is introduced. Finally all facets of confinement are illustrated by a 
simple example of string formation for a quark moving in the field of 
heavy antiquark. 
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1 Introduction 



Ten years ago the author has dehvered lectures on confinement at the XX 
LINP Winter School [1]. Different models of confinement were considered, 
and the idea was stressed that behind the phenomenon of confinement is a 
disorder or stochasticity of vacuum at large distances. 

This phenomenon becomes especially beautiful when one discovers that 
the disorder is due to topologically nontrivial configurations [2] , and the latter 
could be quantum or classical. In that case one should look for appropriate 
classical solutions: instantons, multiinstantons, dyons etc. 

The tedious work of many theorists during last ten years has given an 
additional support for the stochastic mechanism, and has thrown away some 
models popular in the past, e.g. the dielectric vacuum model and the Zi fiux 
model. Those arc reported in [1] (see also refs. cited in [1]). At the same 
time a new and deeper understanding of confinement has grown. 

It was found that the stochastic mechanism reveals itself on the phe- 
nomenological level as a dual Meissner effect, which was suggested as a con- 
finement mechanism 20 years ago [3] . During last years there appeared a lot 
of lattice data in favour of the dual Meissner effect and the so-called method 
of abelian projection was suggested [4], which helps to quantify the analogy 
between the QCD string and the Abrikosov-Nilsen-Olesen (ANO) string. 

At the same time the explicit form of confining configuration - whether 
it is some classical solutions or quantum fiuctuations - is still unclear, some 
candidates are being considered and active work is going on. 

The present review is based on numerous data obtained from lattice cal- 
culations, phenomenology of strong interaction and theoretical studies. 

The structure of the rewiew is the following. After formulation of criteria 
for confining mechanisms in the next chapter, we describe in Section 3 the 
general method of vacuum correlators (MVC) [5] to characterize confinement 
as a property of field correlators. A detailed comparison of superconductivity 
and confinement in the language of MVC and in the effective Lagrangian 
formalism of Ginzburg-Landau type is given in Section 4. 

The Abehan Projection (AP) method is introduced in Section 5 and lat- 
tice data are used to establish the similarity of the QCD vacuum and the 
dual superconductor medium. At the same time the AP method reveals the 
topological properties of confining configurations. 

In Section 6 the most probable candidates for such configurations are 
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looked for among classical solutions. A new principle is introduced to select 
solutions which are able to confine when a dilute gas is formed of those. 

An interconnection of topology and stochastic properties of the QCD 
vacuum is discussed in Section 7. 

In conclusion the search for confining configurations is recapitulated and 
a simple picture of the confining vacuum and string formation is given as 
seen from different points of view. A possible temperature deconfinement 
scenario is shortly discussed. 

2 What is confinement? 

By confinement it is understood the phenomenon of absence in physical 
spectrum those particles (fields) which are present in the fundamental La- 
grangian. In the case of QCD it means that quarks and gluons and in general 
all colored objects cannot exist as separate asymptotic objects. 

To be more exact let us note, that sometimes the massless particles en- 
tering with zero bare masses in the Lagrangian, e.g. photons or gluons, can 
effectively acquire mass - this is the phenomenon of screening in plasma or in 
QCD vacuum above the deconfinement temperature. Due to the definition of 
confinement given above, the screened asymptotic color states cannot exist 
in the confining phase. 

On the other hand in the dcconfined phase quarks and gluons can evolve 
in the Euclidean (or Minkowskian) time sepately, not connected by strings. 
From this fact it follows that the free energy for large temperature t is propor- 
tional to as it usual for the Stefan-Boltzmann gas, and it is natural to call 
this phase the quark-gluon plasma (QGP). It turns out, however, that the 
dynamics of QGP is very far from the ideal gas of quarks and gluons, and a 
new characteristic interaction appears there, called the magnetic confinement 
- this will be discussed in more detail in conclusions, but now we come back 
to the confinement phase. 

It is very inportant to discuss confinement in gauge-invariant terms. Then 
the physical contents of the confinement phenomenon can be best understood 
comparing a gauge-invariant system of an electron and positron (+e, — e) in 
QED and a system of quark antiquark (qq) in QCD when no other charges 
are present. When distance between +e and — e is large the electromag- 
netic interaction becomes negligible and the wave function (w.f.) of {+e, — e) 
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factorizes into a product of individual w.f. (the same is true for a gauge- 
invariant Green's function of (+e, — e). Therefore the notion of isolated e.m. 
charge and its individual dynamics makes sense. In contrast to that in the 
modern picture of confinement in agreement with experiment and lattice 
data (see below), quark and antiquark attract each other with the force of 
approximately 14 ton. Therefore q and q cannot separate and individual 
dynamics, individual w.f. and Green's function of a quark (antiquark) has 
in principle no sense: quark (or antiquark) is confined by its partner. This 
statement can be generalized to include all nonzero color changes, e.g. two 
gluons are confined and never escape each other (later we shall discuss what 
happens when additional color charges come into play). 

It is also clear now that the phenomenon of confinement is connected to 
the formation of a string between color charges - the string gives a constant 
force at large distances; and the dynamics of color charges there without the 
string is inadequate. 

This picture is nicely illustrated by many lattice measurements of poten- 
tial between static quarks. Fig. 1. 

One can see in Fig. 1 clearly the linear growth of potential V{r) — ar at 
large distance r; r > 0.25 fm. 

We shall formulate this as a first property of confinement. 

/. The linear interaction between colored objects. 

To give an exact meaning to this statement, making it possible to check 
in lattice calculations, it is convenient to introduce the so-called Wilson loop 
[6], and to define through it the potential between color changes. This is done 
as follows. Take a heavy quark Q and a heavy antiquark Q and consider a 
process, where the pair QQ is created at some point x, then the pair separates 
at distance r, and after a period of time T the pair annihilates at the point y. 
It is clear that trajectories of Q and Q form a loop C (starting at the point x 
and passing through the point y and finishing at the point x) . The amplitude 
(the Green's function) of such a process according to Quantum Mechanics, is 
proportional to the phase (Schwinger) factor W ~ exp ig J A^jf^d'^x, where 

is the total color vector potential of quarks and vacuum fields, and j/^ - 
the current of the pair QQ, depicting its motion along the loop C. Taking 
into account that is a matrix in color space {A^ — A^T"', a — 1, ...N^ — 1, 
and = ^ct" for SU(2) and T« = for SU(3), where A^^ number of colors 
o"". A" - are Pauli and Gell-Mann matrices respectively) one should insert 
an ordering operator P which orders these matrices along the loop, and the 
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trace operator in color indices tr , since quarks were created and annihilated 
in a white state of the object QQ. Hence one obtains the Wilson operator 
for a given field distribution W{C) — tr P exp ig Jq A^dx^ where we have 
used the point-hke structure of quarks, jij,{x) ~ 5^^\x — x{t))^ 

Now one must take into account, that the vacuum fields A^{x) form a 
stochastic ensemble, and one should average over it. This is a necessary 
consequence of the vector character of Ai^[x) and of the Lorentz-invariance 
of the vacuum, otherwise for any fixed function A^{x) this invariance would 
be violated contrary to the experiments. 

Finally, in field theory in general and dealing with field vacuum in par- 
ticular it is convenient to use the Euclidean space-time for the path-integral 
representation of partition functions, and also for the Monte-Carlo calcu- 
lations on the lattice. There are at least two reasons for that: a technical 
reason; Euclidean path integrals have a real (and positive- definite) measure, 
exp{—SE), where Se is the Eucledian action, and the convergence of the path 
integrals is better founded (but not strictly proved in the continuum). 

Another, and possibly a deeper reason: Till now all nontrivial classical 
solutions in QCD (instantons, dyons etc.) have Eucledian nature, i.e. in the 
usual Minkowskian spacetime they describe some tunneling processes. As 
for the covalent bonds between atoms due to the electron tunneling from one 
atom to another, also the Eucledian configurations in gluodynamics and QCD 
yield attraction, i.e. the lowering of the vacuum energy, and therefore are ad- 
vantageous for the nonperturbative vacuum reconstruction (cf discussion in 
Conclusion). Therefore everywhere in the review we shall use the Eucledian 
space-time, i.e. we go over from Xo,Aq to the real Euclidean components 
X4 = ixo, A4 = UAq. 

Now we come back to the Wilson loop and relate it to the potential. To 
this end we recall that W is an amplitute, or the Green's function of the 
QQ system, and therefore it can be expressed through the Hamiltonian H, 
namely W = exp{~HT), where T is the Eucledian time. For heavy quarks 
the kinetic part of H vanishes and only the QQ potential V{r) is left, so 
that we finally obtain for the averaged (over all vacuum field configurations 
denoted by angular brackets) Wilson loop 



where the loop C can be conveniently chosen as a rectangular r x T, and 




(1) 



5 



confinement corresponds to the linear dependence of the potential V{r), 
V{r) = ar, where a is called the string tension. 

Having said this, one should add necessary elaborations. First, what 
happens when other charges are present, e.g. qq pairs are created from the 
vacuum. In both cases screening occurs, but whereas for ee system nothing 
crucial happens at large distances, the qq system can be split by an additional 
qiqi pair into two neutral systems of qqi and qiq which can now separate - 
the string breaks into two pieces. 

The same is always true for a pair of gluons which can be screened by 
gluon pairs from the vacuum. 

In practice (i.e. in lattice computations) the pair creation from the vac- 
uum is suppressed even for gluons, [7], one of suppression factors for quark 
pairs is 1/Nc [8], another one is numerical and not yet understood, for Nc — 3 
the overall factor is around 0.1, as can be seen e.g. in ratio of resonance 
widths over their masses T/M. 

This circumstance allows one to see on the lattice the almost constant 
force between static q and q up to the distance of around 1 fm or larger. 
This property can be seen in lattice calculations made with the account of 
dynamical fermions (i.e. additional quark pairs), e.g. on Fig.2 one can see 
the persistence of linear confinement in all measiued region with accuracy of 
10%. Therefore we formulate the second property of confinement in QCD, 
which should be obeyed by all realistic theoretical models. 

//. Linear confinement between static quarks persists also in presence of 
qq pairs in the physical region (0.3/m < r < 1.5fm). Another important 
comment concerns quark (or gluon) dynamics at small distances. Perturba- 
tive interaction dominates there because it is singular, which can be seen in 
the one-gluon-exchange force of p-. Comparing this with the confinement 
force mentioned above {a — 0.2GeV'^ = lAton), one can see that perturbative 
dynamics dominates for r < 0.25/m. 

At these distances quarks and color charges in general can be considered 
independent, with essentially perturbative dynamics, which is supported by 
many successes of perturbative QCD. 

Till now only color charges in fundamental representation of SU{Nc) 
(quarks) have been discussed. 

Very surprising results have been obtained for interaction of static charges 
in other representations. For example adjoint charges, which can easily be 
screened by gluons from the vacuum, in lattice calculations are linearly con- 
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fined in the physical region (r < 1.5/m), Fig. 3. 

One can partly understand this property from the point of view of large 
Nf.: the screening part of potential V5 ~ ^exp{—^r) is suppressed by the 
factor 1/Nl as compared with hnear part, but wins at large distances [7] 

< WadAC) >- C^expi-aadjTT) + 9^exp{-V,{r)T) (2) 

c 

The same property holds for other charge representations [10], and more- 
over the string tension (T{j) for a given representation j satisfies an approxi- 
mate relation [7] 



a{fund) C2{fund) 
where C2(j) is the quadratic Casimir operator, 

C2{adj, N,) = Nc, C2ifund, N,) = ~ ^ 



(3) 



2N, 



c 



Correspondingly we formulate the third property of confinement. 

III. Adjoint and other charges are effectively confined in the physical re- 
gion (r < 1.5fm) with string tension satisfying (2). 

We conclude this section with two remarks, one concerning interrelation of 
confinement and gauge invariance, another about many erroneous attempts 
to define confinement through some specific form of nonperturbative gluon 
and quark propagator. 

Firstly, gauge invariance is absolutely necessary to study properly the 
mechanism of confinement. It requires that any gaugc-noninvariant quantity, 
like quark or gluon propagator, 3-gluon vertex etc. vanish when averaged 
over all gauge copies. One can fix the gauge only for the gauge-invariant 
amplitude, otherwise one lacks important part of dynamics, that of confine- 
ment. 

As one popular example one may consider the "nonperturbative" quark 
propagator, which is suggested in a form without poles at real masses. As 
was told above, this propagator has no sense, since on formal level it is 
gauge-noninvariant and vanish upon averaging, and on physical grounds, 
propagator of a colored object cannot be considered separately from other 
colored partner(s), since it is connected by a string to it, and this string 
(confinement) dynamics dominates at large distances. The same can be told 
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about the "nonperturbative" gluon propagator behaving like Xjq^ at small 
g, which in addition imposes wrong singular nonanalytic behaviour of two- 
gluon-glueball Green's function. 

Likewise the well-known Dyson-Schwinger equations (DSE) for one-particle 
Green's functions cannot be used for QCD in the confined phase, because 
again they are i) gaugc-noninvariant ii) there should be a string connected 
to any of propagators of DSE, omitted there. 

More subtle formally, but basically the same is situation with Bethe- 
Salpeter equation (BSE) on the fundamental level, i.e, when the kernel of BSE 
contains colored gluon or quark exchanges: any finite approximation (e.g. 
ladder-type) of the kernel violates gauge invariance and looses confinement. 
When averaging over all vacuum fields is made in the amplitude on the 
other hand, the resulting effective interaction can be treated approximately 
in the framework of BSE. Their value if any might lie in phenomenological 
applications and not in fundamental understanding of the confinement, which 
is the primary purpose of the present review. 

Let us turn now to the order parameters defining the confining phase. To 
distinguish between confined and deconfined phase several order parameters 
are used in absence of dynamical (sea) quarks. One is the introduced above 
Wilson loop (1). Confinement is defined as the phase where the area law 
(linear potential) is valid for large contours, whereas in the deconfining phase 
the perimeter law appears. 

Another, and sometimes practically more convenient for nonzero temper- 
ature is the order parameter called the Polyakov line: 

1 [f^ 
<L{x)>^ - <tr P ex^ig \ ^40^x4 >, I3=l/t (4) 
n Jo 

Since < L > is connected to the free energy F of isolated color charge, 
< L >— exp{—F/T), vanishing of < L > at T < Tc means that F is infinite 
in the confined phase. On the other hand, vanishing of < L > is connected 
to the Z{Nc) symmetry, which is respected in the confining vacuum and 
broken in the deconfined phase. As for the Wilson loops, Polyakov lines can 
be defined both for fundamental and adjoint charges; the first in absence of 
dynamical quarks vanish rigorously in confined phase, while adjoint lines are 
very small there. When dynamical quarks are admitted in the vacuum, both 
Wilson loop and Polyakov line are not order parameters, strictly speaking, 
there is no area law for large enough Wilson loops and < L > does not 
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vanish in the confined phase. However, for large Nc, and practically even 
for A'c = 2, 3 these quantities can be considered as approximate and useful 
order parameters even in presence of dynamical quarks, as will be seen in 
next sections. 

From dynamical point of view the (approximate) validity of area law 
(linear potential) for any color charges even in presence of dynamical quarks 
at distances r < 1.5/m means that strings are formed at these distances and 
string dynamics defines the behaviour of color charges in the most physically 
interesting region. 



3 Field correlators and confinement picture 

As can be seen in the area law (1) the phenomenon of confinement necessarily 
implies appearance of a new mass parameter in the theory - the string tension 
a has dimension [mass]^. 

Since perturbative QCD depends on the mass scale only due to the renor- 
malization through Aqcdj and using to the asymptotic freedom (and rcnor- 
malization group properties) one can express the coupling constant (?(A) at 
the scale A through Aqcd as g'^{A) = , f^o = y-^c ~ |^/) one can 

write 

cr = constAlcD ~ A^exp(- ^^^^^^^ ) (5) 

where A is the cut-off momentum. It is clear from (5) that a cannot be 
obtained from the perturbation series, hence the source of a and of the whole 
confinement phenomenon is purely nonperturbative. Correspondingly one 
should admit in the QCD vacuum the nonperturbative component and split 
the total gluonic vector potential as 

A^^B^ + a,, (6) 

where is nonperturbative and - perturbative part. As for 5^, it can 
be 

a) quasiclassical, i.e. consisting of a superposition of classical solutions like 
instantons, multiinstantons, dyons etc. This possibility will be discussed 
below in Section 6. 
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b) purely quantum (but nonperturbative). The picture of Gaussian stochastic 
vacuum gives an example, which is discussed below. 

It is important to stress at this point, that the formalism of field cor- 
relators, given below in this chapter, is of general character and allows to 
discuss both situations a) and b), quasiclassical and stochastic. In the first 
case, however, some modifications are necessary which will be introduced at 
the end of the chapter. As it was mentioned above, confinement implies the 
string formation between color charges. To understand how string is related 
to field correlators, consider a simple example of a nonrclativistic quark mov- 
ing at a distance r from a heavy antiquark fixed at the origin. As it is known 
from quantum mechanics, the quark Green's function is proportional to the 
phase integral along its trajectory C 

G{f,t) ~< exp ig f A^{r{t'),t')dz^ > 

where the averaging is over all vacuum configurations. It is convenient to 
express though the field strength F^j,, since the latter would be the basic 
stochatic quantities, and this can be done e.g. using the Fock-Schwinger 
gauge 

Jo X 
Hence in the lowest order one obtains 

2 

G{f,t) ~ 1 - y / da^f,{u)da^>^,{u') < F^^{u)F^,^,{u') > +... 

where notations are used such that da^^ = a{u)dx^dui, . 

On the other hand one inroduce the potential V{r) acting on the quark 

G ~ exp{- j V{r, t')dt') ~ 1 - / V{r, t')dt' 

the string formation implies that V{r, t) is proportional to r and this depends, 
as we see, on the field correlators < FF >. 

For the exact Lorcntz-invariant treatment let us introduce gauge^nvariant 
field correlators (FC) and express the Wilson loop average through FC. This 
is done using the nonabelian Stokes theorem [11] 

< W{C) >=< -^Ptr exp ig [ A^dx^, >= (7) 
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<P tr exp ig j^da^,^F^^{u, Zq) > 



Nc 

where we have defined 

rx 

F^,,{u,zo) = ^{zo,u)Fij,^{u)^{u, zo),^{x,y) = Pex.pig / A^dz^ (8) 

and integration in (7) is over the surface S inside the contour C, while Zq is 
an arbitrary point, on which < W{C) > evidently does not depend. In the 
Abelian case the parallel transporters ^{zo, u) and ^{u, zo) cancel and one 
obtains the usual Stokes theorem. 

Note that the nonabelian Stokes theorem, eq. (7), is gauge invariant even 
before averaging over all vacuum configurations - the latter is implied by the 
angular brackets in (7). 

One can now use the cluster expansion theorem [12] to express the r.h.s. 
of (7) in terms of FC, namely [5] 

< W{C) >= — exp V / da{l)da{2)...da{n) < F{l)...F{n) > (9) 

^1 n' J 

where lower indices of da^i, and F^^, are suppressed and F{k) = F^^u^{u'^^\ zq). 

Note an important simplification - the averages <C F{l)...F{n) ^ in the 
color symmetric vacuum are proportional to the unit matrix in color space, 
and the ordinary operator P is not needed any more. 

Eq. (9) expresses Wilson loop in terms of gauge invariant FC, also called 
cumulants [12], defined in terms of FC as follows: 

< F(1)F(2) >=< F(1)F(2) > - < F(l) >< F{2) > (10) 

< F(1)F(2)F(3) >=< F(1)F(2)F(3) > - < F(1)F(2) >< F(3) > - 

- < F(l) X F(2)F(3) > - < F(2) X F(1)F(3) > - < F(l) >< F(2) >< F(3) > 
Let us have a look at the lowest cumulant 

< F{x)<i>{x, Zo)^{zo, y)F{y)^{y, Zo)^{zo, x) > (11) 

It depends not only on x, |/ but also on the arbitrary point zq. In case when a 
classical solution (a dyon or instanton) is present it is convenient to place zq at 
its center, and then zq acquires a clear physical meaning. We shall investigate 
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this case in detail in chapter 6, but now we consider the hmit of stochastic 
vacuum, when the expansion (9) is particularly useful. To this end consider 
parameters on which a generic cumulant <^ F(l)...F(n) ^ depends. When 
all coordinates u'^ coincide with zq, one obtains condensate <C (-P)j„i/,i(0))" 
to which we assign an order of magnitude F". The coordinate dependence 
can be characterized by the gluon correlation length Tg, which is assumed to 
be of the same order of magnitude for all cumulants. Then the series in (9) 
has the following estimate 

< W{C) ^ exp f; M!F'^r2(-i)5 (12) 

where S is the area of the surface inside the contour C. To obtain the result 
(12) we have taken into account that in each cumulant 

« F{x^^\zo)F{x^^\z^)...F{x^^\zo) » (13) 

whenever x and y are close to each other. 

\x -y\<.\x - zq\, \ y- zo \ , 

the dependence on Zq drops out, therefore in (13) in a generic situation where 
all distances | x*^*^ — x'^^'^ |~ Tg <^| a;*^*-' — Zq \ ^ \ x^^^ — Zq \ one can omit 
dependence on zq. 

The expansion in (12) is in powers of (FTg), and when this parameter is 
small, 

FT^ < 1 (14) 

one gets the limit of Gaussian stochastic ensemble where the lowest (quadratic 
in Ffj^iy) cumulant is dominant. 

In the same approximation (e.g. Tg — > while < F^^ > is kept fixed) 
one can neglect in this cumulant the Zq dependence, using the equivalent 
(effective) form of (11) 

D,^uxa = 4r^r < F^^{x)^{x, y)Fxa{y)<^{y, x) > (15) 

The form (15) has a general decomposition in terms of two Lorentz scalar 
functions D{x — y) and Di{x — y) [5] 

Di^uXc = {^nx^ux - -y)+ (16) 

+ 2 ^''■^ [^^^ ■ ^""^ ~ hsCTi^x) + ...]'Di{x - y)} 
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Here the ellipsis implies terms obtained by permutation of indices. It is 
important that the second term on the r.h.s. of (16) is a full derivative by 
construction. 

Insertion of (16) into (9) yields the area law of Wilson loop with the string 
tension a 

< W{C) >= exp{-aS„un) 

j D{x)d''x{l + 0{FT^)) (17) 

where 0{FTg) stands for the contribution of higher cumulants, and Smin is 
the minimal for contour C. 

Note that Di does not enter a, but gives rise to the perimeter term and 
higher order curvature terms. On the other hand the lowest order perturba- 
tive QCD contributes to Di and not to D, namely the one-gluon-exchange 
contribution is 

= ^ (1') 
Nonperturbative parts of D{x) and Di{x) have been computed on the 
lattice [13] using the looking method, which suppresses perturbative fluctu- 
ations, and arc shown in Fig.4. As one can see in Fig.4, both functions are 
well described by an exponent in the measured region, and Di{x) ~ ^D{x) ~ 
exp{—x/Tg) , where Tg ~ 0.2/m. The smallness of Tg as compared to hadron 
size confirms the approximations made before, in particular the stochasticity 
condition (14). One should also take into account that F in (14) is an ef- 
fective field in cumulants, which vanish when vacuum insertion is made and 
therefore can be small as compared with F from the gluonic condensate. 

The representation (16) is valid both for abelian and nonabelian theories, 
and it is interesting whether the area law and nonzero string tension obtained 
in (17) could be valid also for QED (or U{1) in lattice version). To check it 
let us apply the operator ^e^^ap-^ to both sides of (15-16) []. 

In the abelian case, when $ cancel in (15), one obtains (the term with 
Di drops out) 

da < Fap{x)Fxa{y) >= Sxa^f^d^D^x - y) (19) 

where = \ea(3^j.uFfj,u. 

If magnetic monopoles are present in Abelian theory (e.g. Dirac monopoles) 
with the current j^, one has 

daFap{x) = ~jp{x) (20) 
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In absence of magnetic monopoles (for pure QED) the abelian Bianchi iden- 
tity (the second pair of Maxwell equation) requires that 



daFap = (21) 

Thus for QED (without magnetic monopoles) the function D{x) vanishes due 
to (19) and hence confinement is absent, as observed in nature. 

In the lattice version of U{1) magnetic monopoles are present (as lattice 
artefacts) and the lattice formulation of our method would predict the con- 
finement regime with nonzero string tension, as it is observed in Monte-Carlo 
calculations [14]. 

The latter can be connected through D{x) to the correlator of magnetic 
monopole currents. Indeed multiplying both sides of (19) with ^^XajS-^ one 
obtains 

< hi-)jM >= • ^ • ^« - ^4""" " 

The form of Eq.(22) identically satisfies monopole current conservation: ap- 
plying ^ or ^ on both sides of (22) gives zero. 

It is interesting to note, that confinement (nonzero a and D, see Eq. (17)) 
in U{1) theory with monopole currents occurs not due to average monopole 
density < ji{x) >, but rather due to a more subtle feature - the correlator 
of monopole currents (22), which can be nonzero for the configuration where 
< j4{x) >— 0. The latter is fulfilled for the system with equal number of 
monopoles and antimonopoles. 

Note that our analysis here is strictly speaking applicable when stochas- 
ticity condition (14) is fulfilled. Therefore the case of monopoles with Dirac 
quantization condition is out of the region of (14) and needs some elaboration 
to be discussed later in this chapter. 

Let us turn now to the nonabelian case, again assuming stochasticity 
condition (14), so that one can keep only the lowest cumulant (15-16). 

Applying as in the Abelian case the operator l^^ya/Sgf- to the r.h.s. of 
(15-16), one obtains [15] 

< DaF^p{x)^x, y)Fxa{yMy, x) > +^pxa{x, y) = (23) 

= ea/3Xa-Q^D{x -y) 
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The term with Di in (23) drops out as in the Abehan case; the first 
term on the Lh.s. of (23) now contains the nonabelian Bianchi identity 
term which should vanish also in presence of dyons (magnetic monopoles) - 
classical solutions of Yang-Mills theory, i.e. one has 

DaFap{x) = 0. (24) 

It is another question, whether or not in lattice formulation one can vi- 
olate (24) in the definition of lattice artefact monopoles, similarly to the 
Abelian case. We shall discuss this topic when studying lattice results on 
Abehan projected monopoles in Section 5. To conclude discussion of (23) 
one should define A^^^; the latter appears only in the nonabelian case due 
to the shift of the straight-line contour {x.y) of the correlator (15) into the 
position {x + 5x,y), which is implied by differentiation This "contour 
differentiation" is well known in literature [16] and leads to the answer 

A/3Aa(a;, y) = i9 dzpa{z) < tr[Fap{x)(p{x, z)Fap{z)(p{z, y)Fxa{y)(t){y, x)- 
Jy 

(25) 

-Fo,i3{x)(t){x,y)Fxc{y)<t>{y,z)F^p{z)(t){z,x)\ > . 

Especially simple form of (23) occurs when using (23) and tending x to y; 
one obtains [15] 

^|.=o = f < i^:,(0)F,^,(0)F,^jO) > . (26) 

Thus confinement (nonzero a due to nonzero D) occurs in nonabelian case 
due to the purely nonabelian correlator 

< tfFapFp^F^a >— 3 < trEiEjBk > Sijk- 

To see the physical meaning of this correlator, one can visualize magnetic 
and electric field strength lines (FSL) in the space. Each magnetic monopole 
is a source of FSL, whether it is a real object (classical solution or external 
object like Dirac monopole) or lattice artefact. 

In nonabelian theory these lines may form branches, and e.g. electric 
FSL may emit a magnetic FSL at some point, playing the role of magnetic 
monopole at this point. This is what exactly nonzero triple correlator 

< TrFFF > implies. Note that this could be a purely quantum effect, and 
no real magnetic monopoles are necessary for this mechanism of confinement. 
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Till now we have discussed confinement in terms of the lowest cumulant 
- D{x), which is justified when stochasticity condition (14) is fulfilled and 
D{x) gives a dominant contribution. Let us now turn on other terms in 
the cluster expansion (9). It is clear that the general structure of higher 
cumulants is much more complicated than (16), but there is always present 
a Kroneckcr-type term D(a;i, ,X2, .... a;„) J! ^^Mfc similar to D[x — y) in (16) 
and other terms containing derivatives and coordinate differencies like Di. 
The term D{xi,X2, ■■■,Xn) contributes to string tension, and application of 
the same operator ^e^^a3-£- again reveals the nonabelian Bianchi identity 
term (24) and an analog of Apx^ in (23). This means that the string tension 
in the general case is a sum 

oc n « 

a=T a(") , a^") = ^ / < F(l)F(2)...F(n) > da (2)... da (n) (27) 

n=2 ^' 

When the stochasticity condition (14) holds, the lowest term, a^'^\ dominates 
in the sum; in the general case all terms in the sum (27) are important. The 
most important case of such a situation is the discussed above the case of 
quasiclassical vacuum, which we now shortly discuss, shifting the detailed 
discussion to the Section 6. 

For a dilute gas of classical solutions the role of vacuum correlation length 
Tg is played by the size of the solution p. For the correlator < F{xi)...F{xn) > 
the essential nonzero result occurs when all inside the radius p 

of the solution (e.g. dyon or instanton). 

On the other side the typical F of the solution, e.g.: F = F^i,(0), is 
connected to p by the value of topological charge, e.g. for instanton 

Since Q is an integer one immediately obtains that 

{gF-pY--n, n = l,2... (29) 

The same estimate holds for the dyon, solution which is discussed in 
Section 6. 

Hence the series (12) and (27) has a parameter of expansion (gFT^) of 
order of unity and may not converge. A more detailed analysis of the gas of 
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instantons and magnetic monopoles shows that the string tension series (27) 
for instantons (for simphcity centers of instantons and monopoles were taken 
on the plane (12) of the Wilson loop) looks like [17] 

(7 = po(l- < cos/? >); p ^ g j Fi2{z)(fz ^277, (7 = (30) 

while for magnetic monopoles 

/3 = TT, (7 = 2po, (31) 

and po is the surface density of instantons (monopoles). 

Note that a^^^ — po^^^ in both cases is positive, and for instantons 
the total sum for a vanishes while for monopoles cr is nonzero. As we shall 
see below this is in agreement with the statement in chapter 6: all topo- 
logical charges (29) yield flux through Wilson loop equal to 27r(5, and for 
(multi)instantons with Q — n — 1,2, ... Wilson loop is W — exp27r(5i = 1, 
and no confinement results for the dilute gas of such solutions. For magnetic 
monopoles (dyons) topology is different and elementary flux is equal to vr, 
bringing conflnement for the dilute gas of monopoles in agreement with (31). 

Thus the lowest cumulant < F{x)^F{y)^ > might give a misleading 
result in case of quasiclassical vacuum, and one should sum up all the se- 
ries to get the correct answer as in (30-31). Therefore to treat the vacuum 
containing topological charges one should separate the latter and write their 
contribution explicitly, while the rest - quantum fluctuations with FT^ <^ 1 
- can be considered via the lowest cumulants. An example of such vacuum 
with instanton gas and conflning conflgurations was studied in [18] to obtain 
chiral symmetry breaking; this work demonstrates the usefulness of such an 
approach. 

We conclude this chapter with discussion of conflnement for charges in 
higher representations. As it was stated in the previous chapter, our deflni- 
tion of conflnement based on lattice data, requires the linear potential be- 
tween static charges in any representation, with string tension proportional 
to the quadratic Casimir operator. 

Consider therefore the Wilson loop (1) for the charge in some represen- 
tation; the latter was not specifled above in all eqs. leading to (27). One can 
write in general 

A^{x) = AIT", triT^T^) = -5'^'' (32) 
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Similarly to (9) one has for the representation j = (mi,m2...) of the group 
SU{N) with dimension N{j) 

< W{C) = j;^trjexp J2 ^ J da(l)...da(n) < F(l)...F(n) > (33) 

and by the usual arguments one has Eq.(27) . 

Due to the color neutrality of the vacuum each cumulant is proportional 
to the unit matrix in the color space, e.g. for the lowest cumulant one has 

< F(1)F(2) F^(1)F'^(2) > T^X, = (34) 

=< F^(1)F^(2) > j^T:^T^, = M'^Cij) ■ 
where we have used the definition 

T^T^^C2{j)i (35) 
and introduced a constant not depending on representation, 

A(^) ^ ^^J-^ < F%1)F%2) >, (36) 



and also used the color neutrality of the vacuum 



(37) 



For the next - quartic cumulant one has 

< F(1)F(2)F(3)F(4) F"i(l)F"2(2)F''^(3)F'^^(4) > x (38) 

Thus one can see in the quartic cumulant a higher order of quadratic Casimir 
and higher Casimir operators. 

The string tension for the representation j is the coefficient of the diagonal 
element in (34) and (38) 

aU) = C^ij) J ^-^d'x + 0{Cl{3)) (39) 
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where the term 0(C|(j)) contains higher degrees of C2{j) and higher Casimir 
operators. 

Comparing our result (39) with lattice data [10] and Fig.3 one can see that 
the first quadratic cumulant should be dominant as it ensures proportionality 
of a{j) to the quadratic Casimir operator. 

Another interesting and important check of the dominance of bilocal cor- 
relator (of the Gaussian stochasticity) is the calculation of the QCD string 
profile, done in [19]. Here the string profile means the distribution pu of 
the longitudinal component of colorelectric field as a function of distance x± 
to the string axis. This distribution can be expressed through integral of 
functions D{x), Di{x) [18], and with measured values of D, Di from [13] 
one can compute pii(a;j_) and compare it with independent measurements. 
This comparison was done in [19] and shown in Fig. 5. One can see a good 
agreement of the computed pii{x±) (broken line) with "experimental " val- 
ues. Thus MVC gives a good description of data in the simplest (bilocal) 
approximation, even for such delicate characteristics as field distributions 
inside the string. 

4 Dual Meissner mechanism, confinement and 
superconductivity 

The physical essence of the confinement phenomenon is the formation of the 
string between the probing charges introduced into the vacuum, which in 
turn means that the electric field distribution is drastically changes from the 
usual dipole picture (for the empty vacuum) and is focused instead into a 
string picture in the confining vacuum. From the point of view of macroscopic 
electrodynamics of media, this effect can be described introducing dielectric 

— * — * 

function e{x), and electric induction D{x) together with electric field E{x), 
D{x) = e{x)E{x). 

One can then adjust e{x), or better s{D) or e{E) to obtain the string 
formation. Another possibility is to choose the effective action as a function 
of in such a way, as to reproduce string-type distribution. 

This direction was reviewed in [1], the main conclusion reached long ago 
[20] was that the physical vacuum of QCD considered as a medium could be 
called a pure diaelectric, i.e. £ = far from probing charges [21]. 
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It is also shown, that it is possible to choose e{E) in such a way as to 
obtain a string of constant radius [21] or with radius slowly dependent on 
length . We shall not follow these results below, as well as results of the so- 
called dielectric model [1,21], referring the reader to the mentioned literature. 

Instead we focus in this and following chapter on another approach, which 
has proved to be fruitful during last years - confinement as dual Meissner 
effect [3] and ideologically connected to it abelian projection method [4]. 

The physical idea used by 'tHooft and Mandelstam [3] is the analogy 
between the Abrikosov string formation in type II- superconductor between 
magnetic poles and proposed string formation between color electric charges 
in QCD. 

We shall study this analogy here from several different point of view: 

1) energetics of the vacuum - minimal free energy of the vacuum 

2) classical equations of motion -Maxwell and London equations 

3) condensate formation and symmetry breaking 

4) vacuum correlation functions of fields and currents. 

We consider in this chapter the 4d generalization of the Ginzburg-Landau 
model of superconductivity, which is called the Abelian Higgs model with the 
Lagrangian [22] 

^ = -^^i - \DM" - - ^o), D, = d,- teA,. (40) 

This model is known to possess classical solutions-Nielsen-Olesen strings 
, which are 4d generalization of the Abrikosov strings, occurring in type 
II superconductor. The latter can be described by the Ginzburg-Landau 
Lagrangian, when coupling constants are chosen correspondingly [23] . 

The Lagrangian (40) combines two fields: the electromagnetic field {A^, F^,^) 
- which will be analogue of gluonic field of QCD, and the complex Higgs field 
(/9(,t), which describes the amplitude of the Cooper-paired electrons in a su- 
perconductor. When A ^ oo the wave functional \I'{yl^, <^(a;)} has a strong 
maximum around (f{x) = ipo, which means formation of the condensate of 
Cooper pairs of amplitude (po- 

Let us look more closely at the model (40) discussing points l)-4) suc- 
cessively. 

1) The energy density corresponding to the (40) is 

e - + + \Do^\' + ^(l^r - (41) 
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From (41) it is clear that in absence of external sources and for large A the 
lowest (vacuum) state corresponds to ip = ipo = const, and E = B = 0. 
For future comparison with QCD it is worthwhile to stress that formation 
of the condensate < F^^ > is not advantageous for large A, since the mixed 
term l^^v^P will give very large positive contribution: condensate of electric 
charges (/? repels and suppresses everywhere electromagnetic field F^^^. The 
same situation occurs in the nonabelian Higgs - the Georgi - Glashow model: 
there appear two phases depending on values of A, g and for large A the 
deconfined phase with (p — (po persists. Here the word deconfined means that 
external (color) electric charges are not confined, while magnetic monopoles 
can be confined. For smaller values of A one can reach a region where it is 
advantageous to form the condensate < F^^ > and then possibly also nonzero 
D{x). As we discussed in Section 3 then there is a possible confinement of 
colorelectric charges in this phase. But let us come back to the Abelian Higgs 
model and the limit of large A, which is of primary interest to us. 

In this case it is advantageous to form the condensate of electric charges 
{(p{x) = (fio), condensate of electromagnetic field is suppressed, and one ob- 
tains Abrikosov- Nielsen-Olesen (ANO) strings connecting magnetic poles - 
this is the confinement of magnetic monopoles and mass generation of e.m. 
field - leading to the deconfinement of electric charges. 

In the dual picture one would assume that condensate of magnetic charges 
(monopoles) would help to create strings of electric field, connecting (color) 
electric charges, yielding confinement of the latter. Our tentative guess is that 
the condensate of magnetic monopoles (dyons) in gluodynamics is associated 
with the gluonic condensate < F^^^ >. 

2) We now discuss the structure of the ANO strings from the point of 
view of classical equations of motion. The Maxwell equations are 



where j is the microscopic current of electric charges, including the conden- 
sate of Cooper pairs. 

To obtain a closed equation for B one needs the specific feature of super- 
conductor in the form of Londons equation 



The latter can be derived from the Ginzburg-Landau-type Lagrangian 
(40). Indeed writing the current for the Lagrangian (40) in the usual way. 



— * — * 

rotB — j 



(42) 



rotj = S-^B 



(43) 
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we have 

J, = te{^+d,^ - - 2eM^|(/)p (44) 

Assume that there exists a domain where 4> is aheady constant, and is 
still nonzero (we shall define this region later in a better way) , and apply the 
rot operation to both sides of (44). Then one obtains the London equations 
(43) and S is defined by 

S-^ = 2e^(f>l (45) 
Insertion of (43) into (42) yields equation for B 

AB - S-^B = (46) 

— * 

Solving (46) one obtains the exponential fall-off of B away from the center 
of the ANO string, 

B{r) = constKo{r/5), B{r) ~ exp{-r/5), r > 5 (47) 

From (47) it is clear that is thephoton mass, generated by the Higgs 
mechanism. On the other hand the field has its own correlation length ^, 
connected to the mass of quanta of the field (the " Higgs mass" ) 

e = l/m^, ml ^2X4 (48) 

As it is known [23] the London's limit for the superconductor of the second 
kind, corresponds to the relations 

5 > e < A (49) 

One can also calculate the string energy per unit length-thc string tension - 
for the string of the minimal magnetic fiux 27r. Using (41) and (47-48) one 
obtains [23] 

7r S 

(TANO^^ln-, VC>1 (50) 

Note that the main contribution to aANO (50) comes from the term |-D0p in 
(41). 

Thus the physical picture of the ANO string in htc London's limit implies 
the mass generation of the magnetic quanta, m — 1/6, which are much less 
than the Higgs mass m<^ = 1/^- 
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It is instructive to see how the screening of the magnetic field (mass 
generation) occurs: 

First magnetic field creates around its fiux (field-strength fine) the circle 
of the current j from the superconducting medium, as described by London's 
equations (43). Then Maxwell equation (42) tells that around the induced 
current j there appears a circulating magnetic field, directed opposite to the 

original magnetic field B and proportional to it. As a result magnetic field 
is partly screened in the middle and finally completely screened far from the 
center of the magnetic fiux (which is the ANO string) . 

The profile of the string, i.e. i?ii(r) as a function of the distance of the 
string axis is shown in Fig.6, as obtained from Eqs. (47). 

It is interesting to compare this ANO~string profile with the correspond- 
ing profile obtained in the gluodynamics. As was discussed in the previous 
chapter, the distribution of the parallel component of the colorelectric field 
was measured in [24] and is also exponentially decreasing far off axis, as can 
be seen in Fig.6, and both profiles of the ANO string and of the QCD string 
in Fig.6 are very much similar, thus supporting the idea of dual Meissner 
mechanism . 

To look more closely at the similarity of classical equations (4.9-4.10) 

with the corresponding equations, obtained in lattice Monte-Carlo simula- 
tions, one needs first an instrument to recognize the similarity of effective 
Lagrangians in the model (40) and in QCD. This is discussed in the next 
chapter. 

We conclude point 2) of this chapter discussing parameters Sqcd and 
Cqcd, parameters dual-analogical to S and { of (44-45). This is done via the 
Abelian projection method in [25] with the result (see Fig. 7) 

Sqcd ~ ^qcd ~ 0.2/m (51) 

Thus in QCD situation is somewhat in the middle between the type I and 
type II dual superconductor. The calculation of the effective potential V{(j)) 
in SU(2) gluodynamics, made recently in [26], shows a two-well structure, 
but with a rather shallow well, hence the effective Xqcd in the Lagrangian 
(40) is not large, again in agreement with (51). A very interesting discussion 
of properties of dual monopoles and their measurement on the lattice is 
contained in [27]. 

3) Condensate formation and symmetry breaking 
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The phenomenon of superconductivity is usually associated with the for- 
mation of the condensate of Cooper pairs (although it is not necessary). 

The notion of condensate is most clearly understood for the noninteract- 
ing Bose-Einstein gas at almost zero temperature, where the phenomenon of 
the Bose- Einstein condensation takes place. Ideally in quantum-mechanical 
systems condensate can be considered as a coherent state . 

When interaction is taken into account, the meaning of the condensate is 
less clear [28]. 

In quantum field theory one associates condensate with the properties of 
the wave functional and/or Fock columns. Again in the case of no interaction 
one can construct the coherent state in the second-quantized formalism (see 
e.g. the definition of wave operators in the superfluidity case in [29]) 

One of the properties of such a state is the fixed phase of the wave func- 
tional, which means that the U{1) symmetry is violated. 

The simplest example is given by the Ginzburg-Landau theory (40), 
where in the approximation when A oo, the wave functional ^{<^} can 
be approximated by the classical solution if = ipo, ^'{9?} — > "iflipo}- 

The solution (p — (po, where </?o has a fixed phase, violates U(l) sym- 
metry of the Lagrangian (40), and one has the phenomenon of spontaneous 
symmetry breaking (SSB) [30]. The easiest way to exemplify the SSB is the 
double-well Higgs potential like in (40). Therefore if one looks for the dual 
Meissner mechanism in QCD (gluodynamics) one may identify the magnetic 
monopole condensate (f, dual to the Cooper pair condensate (p — (fo and find 
the effective potential V{(f), demonstrating that it has a typical double-well 
shape. Such an analysis was performed in [26] using the Abelian projection 
method and will be discussed in the next chapter. 

4) Finally in this chapter we shall look at the dual Meissner mechanism 
from the point of view of field and current correlators [31]. This will enable 
us to formulate the mechanism in the most general form, valid both in (quasi) 
classical and quantum vacuum. 

One can use the correlator (16) to study confinement of both magnetic 
and electric charges (the latter was discussed in chapter 3 - Eq. (19) and 
subsequent text). Let us rewrite (16) for correlators of electric and magnetic 
fields separately 

f)r)E fifjE 
< E,{x)E,{y) >= 5,,{D^ + + hl—^) + h,h,-^ (52) 
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< H,{x)Hj{y) >= 5,j{D^ + D^ + h'^) - h,hj-^ (53) 

where h^ = x^,- y^, = h^h^. 

In (52)- (53) we have specified correlators Di for electric and magnetic 
fields separately, since in Lorentz-invariant vacuum = , Df — , 
otherwise, e.g. in the Ginzburg-Landau model, electric and magnetic corre- 
lators may differ, as also in any theory for nonzero temperature . 

Now let us compare the Wilson loop averages for electric and magnetic 
charges. In case of electric charges the result is the area law (17) with the 
function D —>■ responsible for confinement. 

Now consider a magnetic charge in the contour C in the 14 plane, the 
corresponding Wilson loop is 

< W{C) >=< exp{ig j Fudau) >= exp{-a*Smin) (54) 

Here F^,, is the dual field, F^,, = \e^,vai3Fai3, and = Hi. Prom (17) one 
obtains 

<7* = ^ / d'xD^{x){l + 0{T'^H)) (55) 

Eq. (55) is kind of surprise. For electric charges yields confinement 
and is nonperturbative, supported by magnetic monopoles, cf. Eq. (22), 
while contributes perimeter correction to the area law and contains also 
perturbative contributions like Coulomb term. 

The same electric charges in the (23) - plane Wilson loop bring about 
again the area law (17) with D — > . Also is coupled to magnetic 
monopoles in the vacuum, indeed taking divergence from both sides of (53) 
one gets 

< divH{x), divH{y) >= -d^G^{x - y) (56) 

Instead confinement of external magnetic charges, Eq. (4.16), is coupled to 
the function D'^ (or Df if one takes < W{C) > for the loop in the (23) 
-plane) . 

Thus duality of electric - magnetic extemal charges requires interchange 

At this point one should be careful and separate out perturbative inter- 
action which is contained in Di, namely one should replace in (55) by 
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where 

= - ^ (57) 

It is important to stress again that it is the nonperturbative contents of 
correlators which may create a new mass parameter hke a* and which should 
enter therefore in correlators (the perturbative term contributes the usual 
Coulomb-like interaction which is technically easier consider not as a part of 
, but to separate at earlier stage, see e.g. [32]). 

This mass creation can be visualized in the Ginzburg-Landau model, 
where Di can be computed explicitly from (40) 

D^^ix -y) = (e^l^p - d%' ^ e-l^-^l (58) 

where in the asymptotic region = \(J)q\ and m = e|0|o = 1/S. 

From the point of view of duality D^{x) Eq.(49) should be compared 
with the behaviour of D^{x), which was discussed above in Section 3 and 
measured in [13]. 

D{x) = D^{x) = e-'^l^l, f« IGeV (59) 

Eqs. (58-59) demonstrate validity of dual Meissner mechanism on the level 
of field correlators. 



5 Abelian projection method 

In the previous section it was shown how the string is formed between mag- 
netic sources in the vacuum described by the Abelian Higgs model. 

In the dual Meissner mechanism of confinement [3] it is assumed that the 
condensate of magnetic monopoles occurs in QCD, which creates strings be- 
tween color electric charges. There are two possible ways to proceed from this 
point. In the first one may assume the form of Abelian or nonabelian Higgs 
model for dual gluonic field and scalar field of magnetic monopoles. This 
type of approach was pursued in [33] and phenomenologically is quite suc- 
cessful: the linear confinement appears naturally and even spin-dependent 
forces are predicted in reasonable agreement with experiment [34]. We shall 
not go into details of this very interesting approach, referring the reader to 
the cited papers, since here lacks the most fundamental part of the problem 
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- the derivation of this dual Meissner model from the first principles - the 
QCD Lagrangian. 

Instead we turn to another direction which was pursued very intensively 
the last 8-10 years - the Abelian projection method dating from the seminal 
paper by 'tHooft [4]. 

The main problem -how to recognize configurations with properties of 
magnetic monopoles, which are responsible for confinement. The 'tHooft's 
suggestion [4] is to chose a specific gauge, where monopole degrees of free- 
dom hidden in a given configuration become evident. The corresponding 
procedure was elaborated both in continuum and on lattice [35] and most 
subsequent efforts have been devoted to practical separation (abelian projec- 
tion) of lattice configurations and study of the separated degrees of freedom, 
and construction of the effective Lagrangian for them. We start with the 
formal procedure in continuum for SU(N) gluodinamics, following [4,35]. 

For any composite field X transforming as an adjoint representation, like 

X ^X' = VXV-^ (60) 
let us find the specific unitary matrix V (the gauge), where X is diagonal 

X' = VXV-'^ = diag{\i, A2, ...A^v) (61) 

For X from the Lie algebra of SU(N), one can choose Ai < A2 < A3 < 
...Ajv- It is clear that V is determined up to the left multiplication by a 
diagonal SU(N) matrix. 

This matrix belongs to the Cartan or largest Abelian subgroup of SU(N), 
V(l)^-i C SU{N). 

Now we transform to the gauge (61) 

A^^V{A^+'-d;)V-^ (62) 

and consider how components of A^ transform under U{1)^~^ . The diagonal 
ones 

= (63) 

transform as "photons": 

< ^ < = < + (64) 
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while nondiagonal, = (^4^, transform as charged fields. 

= exv{i{a, - a,)]Cj- (65) 

But as 'tHooft remarks [4], this is not the whole story - there appear 
singularities due to a possible coincidence of two or more eigenvalues Aj, and 
those bear properties of magnetic monopoles. To make it explicit consider 
as in [35] the "photon" field strength 

^VF^^V-^ ^ig^{A^^-d;)V-\ V{A,^-d,)V-^\ (66) 
and define the monopole current 

Since F^^^ is regular, the only singularity giving rise to i^T^ is the commutator 
term in (66), otherwise the smooth part of does not contribute to 
because of the antisymmetric tensor. 

Hence one can define the magnetic charge m'(Q) in the 3d region Q, 



m 



(Q) = / = ^ / ciV^.£^.,,/;, (68) 

JO. OTT JdVL 



Consider now the situation when two eigenvalues of (61) coincide, e.g. 
Ai = A2. This may happem at one 3d point in Q , i.e. on the line in 
the 4d, which one can visualize as the magnetic monopole world line. The 
contribution to m*(f2) comes only from the infinitesimal neighborhood of 

mX5e(x(^^)) = / d\^,E^,p,^dpV-\Vd,V-% = 

= / d^o^,e^,,,dp^d,V%, (69) 

The term Vd„V~^ is singular and should be treated with care. To make it 
explicit, one can write. 

y ^w{^ co^i^ + i^^4>sin\, ^ ^^^^ 
28 



where is a smooth SU(N) function near x^^\ Inserting it in (69) one 
obtains 

m'(S,(a:«)) = ^ / d^a^,e^,p,dp{l-cos9)d,cl>[a3]u (71) 

where and 6 are azimuthal and polar anges. 

The integrand in (71) is a Jacobian displaying a mapping from S^{x^^^) 
to (9,(1)) ^ SU{2)/U{1). 

Since 

n2{SU{2)/U{l)) = Z (72) 

the magnetic charge is m* = 0, ±1/2, ±1, ... 

From the derivation above it is clear, that the point x = x^^\ where 
Ai(,x'^^)) = A2(.t'-^^), is a singular point of the gauge transformed and a*^, 
and the latter behaves near x = x^^^ as 0( ); and abelian projected field 

strength /^^ is 0( |^_^(i)|2 ); like the field of a point-like magnetic monopole. 

However several points should be stressed now 

i) the original vector potential and F^^, are smooth and do not show 

any singular behaviour 

ii) at large distances /^j, is not, generally speaking, monopole-like, i.e. 
does not decrease as |^_^(i)|2 ) So that similarity to the magnetic monopole 
can be seen only in topological properties in the vicitity of the singular point 
xW. 

iii) the fields A^, have in general nothing to do with classical solutions, 
and may be quantum fluctuations. Actually almost any flcld distribution in 
the vacuum may be abelian projected into a^, /^^ and magnetic monopoles 
then can be detected. 

Examples of this statement will be given below, but before doing that 
we must say several words about the choice of the field X in (60) and more 
generally about the explicit gauge choice. 

By now the most popular choices for the adjoint operator X, (60) is for 
nonzero temperature the fundamental Polyakov line 

Lij{x) — {Pexpig dx4A4{x4,x))ij (73) 

%f 

where P = 1/T, and T is the temperature, A^{x4^, x) is required to be periodic 
in X4 and i,j are fundamental color indices. 
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Another widely used gauge is the so-called maximal abelian gauge (MAG) 
[35] , which in lattice notations can be expressed as a gauge where the follow- 
ing quantity is maximized 

R = Y.Tr{asU,{s)a,U+{s)) (74) 

Here the link matrix ~ expigA^^^Azn, and U is gauge transformed with 
the help of V 

U^{s)=Vis)U^{s)V-\s + fi) (75) 

In the continuum MAG is characterized by the condition, which in the SU(2) 
case looks most simply 

(d.TigADA^^O, A^^Al±zAl (76) 

As will be seen the choice of gauge in the abelian projection method is crucial: 
e.g. for the minimal abelian gauge, corresponding to the minimum of it! (74), 
Abelian projected monopoles have no influence on confinement [36]. 

Since the total of papers on abelian projection is now enormous, let us 
discuss shortly the main ideas and results. The most part of results are 
obtained doing abelian projection on the lattice. (A short introduction and 
discussion of lattice technic is given in [27]). The separation of monopole 
degrees of freedom was done as follows. For an abelian projected link (75) 
one can define the U{1) angle 9^ 

Un = expiiO^a-i), -tt <9^<7r (77) 

Then for the plaquctte f/^j, ~ U^Uj^Y^ exp^iO^i^a^) one can write 4n < 

9ij,u < 47r and define the "Coulomb part" of Ofj,^, i9^i, = mod2Tr{d^v}-, so that 

V = ^^^^ + 2?™^, (78) 

and rijuj, counts the number of Dirac strings across the plaquctte (/xz^). 

Now one can calculate different observables in AP and find the contribu- 
tuion to them separately of the " Coulomb part" Ofj,^ and the monopole part 
m^iy = 27171 This was done in a series of papers and the monopole dom- 
inance was demonstrated i) for the string tension [38] ii) for fermion prop- 
agators and hadron masses [39] iii) for the topological susceptibihty [40] . 
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In Fig. 8 a comparison is made of the total AP contribution to the string 
tension and of the monopole part, which is seen to dominate as compared to 
the Coulomb part. 

Another line of activity in AP is the derivation of effective Lagrangians 
for AP degrees of freedom. ( See e.g. in [41]). The resulting Lagrangian how- 
ever docs not confine AP neutral objects like "photons" and the latter should 
contribute to the spectrum in contradiction to the experiment. Therefore the 
Lagrangian is not very useful both phenomenologically and fundamentally. 
We shall not discuss these Lagrangians referring the reader to the cited liter- 
ature. Let us come back to the investigation of the confinement mechanism 
with the help of the AP method. 

A direct check of the dual Meissner mechanism of confinement should 
contain at least two elements: a detection of dual London current and a 
check of magnetic monopole condensation. The first was done in [42]. The 
dual of the London equation (43) is 

E^5\otjM, S = l/m (79) 

For the ideal type -II superconducting picture one needs that 6 ^ ^, where 
^ = m^^ is the Higgs mass (corresponding to the magnetic monopole con- 
densation). In practice in [42] it was found that 5 ~ ^-iid therefore the 
condensate is " soft" , implying that exact solutions of ANO string are neces- 
sary to compare with. Such analysis was done in [43] for SU{2) and SU{3) 
and recently in [44]. 

The string profile p{x±) was also studied using AP, and the resulting 
p{x±) [45] is similar to the one, obtained in the full lattice simulation [46]. 
Thus the whole picture of currents and density is compatible with the dual 
superconductivity. 

Now we come to the second check - of monopole condensation. Of special 
interest is the problem of definition of monopoles. In the first lattice studies 
[35] the AP monopoles in the maximally Abelian gauge have been indentified 
through magnetic currents on the links of the dual lattice, and the perimeter 
density of the currents was measured below and above transition temperature 
Tc, showing a strong decrease of this density at T > T^. 

Later it was realized [47] that monopole density defined in this way may 
not be a good characteristics of dual superconductivity and monopole con- 
densate, and for the latter one needs the creation operator of the magnetic 
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monopolc. In this way one can define the dual analog of the Higgs field 
and settle the question of condensation. The general mathematical construc- 
tion of the monopole creation operator was given before in [48], and several 
papers [49] used this U(l) construction for the AP monopoles. 

Another construction was used for the U(l) theory in [50] and as in [49] 
the condensation of monopoles was also demonstrated. In the SU(2) case the 
analysis was performed in [51] and [26]. 

In the last paper the important step was done in finding an effective 
potential V{(p) for the monopole cration operator ip (defined similarly to 
Froelich-Marchetti [48], the exact equivalence to [48] is still not proved in 
[26]). 

If one believes in the dual Meissner picture of confinement, then one 
should expect the two-well structure of V{(p), symmetric with respect to the 
change if —ip. One can sec in Fig.9a the r.h.s. of V{(j)) (for positive 
0) which indeed has a minimum at = (f)c, shifted to the right from = 
0, for values of /5 = p in the region of confinement. In Fig. 9b the same 
quantity ^(0) is measured in the deconfinement phase, and as one can see, 
the minimum ip = ipc is at zero, 0c = 0. In this way the analysis of [26] gives 
an evidence of AP monopole condensation. Note, however, that strictly 
speaking the condensation should be proved in the London limit (A oo), 
otherwise quantum fluctuations of (/? might prevail for a shallow well like 
that in Fig. 9a (we remind that the system has a finite number of degrees of 
freedom in a finite lattice volume) . Till now we have said nothing about the 
nature of configurations which due to AP disclose the magnetic monopole 
structure and ensure confinement. They could be classical configurations 
or quantum fluctuations (the last possibility is preferred by the majority of 
researches) . 

Recently an interesting AP analysis of classical configurations has been 
done [52-54]. We shall mostly pay attention to the first paper [52], created 
a wave of further activity. The authors of [52] make AP analytically of an 
isolated instanton, multiinstanton and Prasad-Sommerfield monopole and 
demonstrate in all cases the appearance of a straight-line monopole current. 
In the first case the current is concentrated at the centre of instanton and its 
direction depends on the parametrization chosen. Later on in [53] the same 
type of analysis was made numerically with the instanton-antiinstanton gas, 
and the appearence of monopole-current loops was demonstrated, of the size 
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of instanton radius and stable with respect to quantum fluctuations. 

Hence everything looks as if there is a hidden magnetic monopole inside 
an instanton. This result is extremely surprising from several points of view. 
First of all the flux of magnetic monopole through the Wilson loop is equal to 
TT (this will be shown in the next Section), whereas the same flux of instanton 
is equal to zero (modulo 27r), this fact helps to understand why the monopole 
gas may ensures conflnemet, while the instanton gas does not. Therefore the 
identification of monopoles and instantons is not possible. Secondly, confine- 
ment in the instanton gas was shown to be absent in several independent 
calculations [55], and the appearance of rather large monopole loops in this 
gas [53] looks suspiciously. To understand what happens in the AP method, 
and whether it can generate monopoles where they were originally absent, 
we turn again to the case of one instanton [52] and take into account that 
AP contains a singular gauge transformation, which transforms an originally 
smooth field F^,^ into a singular one, namely [52] 

F^,{x) = VF,,{x)V-' + (80) 

where 

F-"^(x) = '-V{x)[d^d, - d,d^]V-\x) (81) 

The matrix V{x) is singular and F^^^ therefore does not vanish. Hence the 
correlator function D{x), of the fields (15) and currents (22) aquires due to 
this singular gauge transformation a new term 

D(x) D{x) + D"*"^(a;), (82) 

where 

D''''s{x -y)^< F^*"^(x)F*'"^(7/) > (83) 

The same F^"'^ causes the appearence of a magnetic monopole current which 
passes through the center of an isolated instanton [52] . Therefore the AP 
method indeed inserts a singularity of magnetic monopole type into practi- 
cally any configuration and therefore cannot be a reliable method of separa- 
tion of confining configurations. On the other side the use of AP for lattice 
configurations reproduces well the conflnement observables, e.g. the string 
tension [38], which means that confining contributions pass the AP test very 
well. 
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From this point of view it is interesting to consider another example - 
the dyonic (Prasad-Sommerfield) solution, also studied in [52]. It appears 
that the AP monopole current passes exactly through the multiinstanton 
centers and coincides with the trajectory of the physical dyon, which can be 
calculated independently, i.e. the AP magnetic monopole current coincides 
with the total magnetic monopole current. This again supports the idea that 
the genuine confining configurations are well projected by the AP method. 

6 Classical solutions as a possible source of 
confinement 

In the last Section we have seen, that the AP method gives no answer to the 
question what are the confining configurations, they could be both classical 
fields and quantum fluctuations. Nothing about it can be said till now from 
the field correlators, however lattice measurements of correlators yield some 
information on the possible profile of confining configurations. 

At the same time an intersting information can be obtained from the 
lattice calculations using the so-called cooling method [56], where at each 
step of cooling quantum fluctuations are suppressed more and more, and 
configurations evolve in the direction of decreasing of the action. Roughly 
these results demonstrate that out of tens of thousands of original configu- 
rations (which are mostly quantum noise) at some step of cooling only few 
(sometimes 15-25) are retained, which ensure the same string tension, as the 
original - "hot " vacuum. With more cooling the mimbcr of configurations 
drops to several units (they are mostly (anti)instantons) and confinement 
disappears. 

Thus one can think that confining configurations in the vacuum differ 
from usual quantum fluctuations, and their action is probably larger than 

instantonic action or they arc less stable. 

Therefore it is interesting to look carefully into all existing classical solu- 
tions and check whether one can build up the confining vacuum out of those 
solutions. 

In this chapter we shall study several classical solutions: instantons, dyons 
and lattice periodic instantons, and give a short discussion of some other 
objects, like torons. 
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After disposing the individual properties of those we specifically concen- 
trate on the contribution of each of the object to the Wilson loop- what we 
call the elementary flux of the object - and argue that flux proportional to 
TT of dyons and twisted instantons with Q = 1/2 may yield confinement, in 
contrast to the case of instantons with fiux equal to 27r. 

To prove this one should construct the dilute gas of objects, which we do 
in the most nontrivial example of dyons. 



6.1. Classical solutions i.e. solutions of the equation 

D^F^, = (84) 

can be written in the ( ant i) self dual case in the form of the so-called 'tHooft 
ansatz [57] or in the most general form of Atiya-Drinfeld-Hitchin-Manin 
(ADHM) [58]. 

In the first (simpler and less general) case one has 

K = ~%.dMW, (85) 

where 77^^^ is the 'tHooft symbol, 

%v = (^a^Lu^^i, V = 1,2, ?> and 5a^.,/x = 4 or - 5^^, z/ = 4, (86) 
and W satisfies equation ff^W — with a particular solution 

N 2 

W^l + Y] (87) 

Here p„, = 1, ..N are real parameters. In the simplest case, = 1, one 
has the instanton solution [59] with size p„ = p and position a;^^-* . For p„ finite 
and A^ arbitrary (85), (87) gives a multiinstanton solution with topological 
charge Q = N. In particular for A^ — > 00, x^^^ = nb, x^-^^ = r^, one gets 
the Harrington-Shepard periodic instanton [60], for which periodically 
depends on (Euclidean) time 0:4 = t. 

A specific feature of instantons is their finite size: fields F^^ fall off at 
large distances from the center as a;"'*. This is very different from the case 
of magnetic monopole, where fields decay as a;~^. 
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Another class of solutions with these latter properties can be obtained 
from (85), (87) in the limit p„ = p — oo. 

The case oi N — 1, when W — ^^_^(o))2 > yields no solution, since it is a 
pure gauge. 

The next case, = 2, is gauge equivalent to the (anti)instanton with the 
position ^{x^^^ + x^'^^) and size ^\x^^^ — x'^'^^\. 

We shall be interested in the case when — > oo,pi = p— >• oo, x^^^ — 
nb, x^"^ = R and call this solution a dyon, since, as we demonstrate below, it 
has both electric and magnetic long-distance field. 

In this case W can be written as 

oo 1 

W= J2 = ' (88) 

n=-oo {x- Ry + {x4 - nby 

and one can rewrite (88) using variables 

7|f - R\=r, X47 = t, 7 = 2n/b (89) 

as 

W^^-^^. (90) 
2r car — cost 

Prom (85) one obtains vector-potentials 

A A 1 ^hr , 5iaSint , 

Aia = BaikTiki- - cthr + — -) -, (91) 

r car — cost car — cost 

1 shv 
Aia = na{- - Cthr + -), (92) 

r cnr — cost 

with 

{x- R)/\x- R\. 

One can notice, that the dyonic field in this (singular, or 'tHooft) gauge 
is now long ranged in spatial coordinates 

A,a-\, F,.^^ (93) 

and periodic in "time" t. 
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Even more similarity to the magnetic monopole field can be seen, when 
one makes the (singular) gauge tranformation [61] 

i^ = C/+(A^ + ^9^)C/, C/ = exp(^^), (94) 



with 

tgO = Wi 



r 



In this gauge, sometimes called Rossi gauge, looks hke Prasad-Sommerfield 
solution [62]: 

Aa^ f{r)en,anb, /(O = — (1 - ^), (95) 

gr shr 

Ma = v{r)na, <p{r) = —{rcthr - 1). (96) 

gr 

Note that does not depend on time, it describes a static dyonic solu- 
tion, since it has both (color) electric and (color) magnetic field: 

Eka = Bka = Saki-f ' f M + nafikif - f/v + gf). (97) 

One may further gauge rotate Eka-i Bka to the quasiabelian gauge, where the 
only long-ranged component is along 3 axis: 

E'ks = B^r ^ oo) ~ -^n,. (98) 

Eq.(98) justifies our use of the name dyon for the solution and demonstrates 

its clear similarity to the magnetic monopole. Note also, that A^^^ (96) is 
tending to a constant at spacial infinity, like the Higgs field component of 
the 'tHooft-Polyakov monopole [63]. 

The total action of the dyon is calculated from (95), (96) or (91), (92) to 

be 

S-\j d'r£ dt{Bl, + El,) = ^T, (99) 

where T is the length of the " dyonic string" , in terms of the number of 
centers in (87) it is T = b{N — 1). For the given N one also has 

S{N)^^Q{N), Q{N)^N-1. (100) 
9 
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6.2. This short section is devoted to another type of classical solutions - 
those depending on boundary conditions and defined in finite volume. Here 
we consider torons and instantons on torus [64], which obey the twisted 
boundary conditions (b.c.) in the box < < a^. Periodic b.c. are 
imposed modulo gauge transformation (twisted b.c.) 



^x{x^, = a^) = 



A,{x, = 0)-z-^ (101) 



To ensure selfconsistency of ^4^ on the lines, four functions n^(/x = 1, ...,4) 
should satisfy conditions 

ni{x2 = 02)^^2(2:1 = 0) = 1^2(2:1 = ai)ni{x2 = 0)Zi2, (102) 

and analogous conditions for 1, 2 — > where Zi2eZ{N) is the center of the 
group SU(N), 

Zi^u = exp(2ni^^ , n^^ = -n^^. (103) 

Here n^^ are integers not depending on coordinates x^. 

The twisted solutions ^4^ (101) contribute to the topological charge 



f Tr{F^,F^,)d^x = - ^, (104) 

J\x,A<a,, ^ >^ ' TV 



IGtT^ J\x^\<a^ 



where - an integer and x = In^^n^^ = ni2nu + ?T'13?T'42 + ?T'14?T'23- 
The action in the box is bounded from below 



ljTr{F,^F,^)d'x>^\i.-^\. (105) 



Consider e.g. the case of 7134 = —ni2 — 1, all other Uf^i, = 0, x = 1. Then 
one has 

^a(^) = --E^^> = (106) 
where = ^m] = k + l = N, 

uj — 2t: • diag{l, I, —k, —k), (107) 
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the matrix oo has k elements equal to / and / elements equal to (—A;), and 
a condition is imposed 0102(0304)"^ ~ i ~ ^ simple example take 

SU(2) and cubic box, then k — I — l,u; — 27rT3 and 

Ax{x) = — -— af,xx^, ai2 = 0:34 = 1- (108) 

This solution is selfdual and the following relation holds 

W,„ = r.F,„F„„ = j^-^. (109) 

Prom (109) one can see that toron (108) is a particular case of a selfdual 
solution with constant field F^jy 

A^{x) = F^^xJ^, (110) 

where the amplitude of constant field F^^, is quantized. For constant (anti)self- 
dual field the analysis of Leutwyler [65] tells that such solutions are stable 
with respect to quantum fluctuations. 

The flux through the Wilson loop for the solution (108) is the planes (12) 
or (34) is 

f S 
Pexp{ig / Af,dx,j,) = exp{-i'K —t^) , (111) 
J c a 

where S is the area, bounded by the contour C. 

As we shall see in the next section the flux equal to vr for S = \s a, 
property very important for confinement. Another interesting property of 
torons, not shared by any other solutions, that its action is proportional to 
I/Ncq"^ and therefore stays constant for large Nc, where — Qq/Nc- We 
come back to this property in conclusions. 

Another type of twisted solutions arc twisted instantons [64]. These are 
solutions with topological charge Q (104) and u integer non zero. 

Those solutions have been seen on the lattice [66] , and the profile (distri- 
bution of trF^^{x)) is very close to that of the usual instanton. 

Unfortunately the analytic form of twisted instanton is still unknown; the 
top charge was found to be 1/2 [66], and the extrapolated string tension is 
probably nonzero. These two facts are not accidental - in the next section 
we show that halfinteger top. charge ensures a flux of n and that in turn 
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may lead to confinement. 



6.3. In this subsection we compute elementary flux inside Wilson loop 
for (multi)instanton, dyon and twisted instanton and connect properties of 
elementary flux to confinement in the gas of classical solutions [67]. 

Consider a circular Wilson loop in the plane (12) and take in the form 
of the 'tHooft ansatz (5.2), where is fixed and xf'-* = 0,x^^^ = nb. In 
this way one can study the case of instanton {N = 1), periodic Harrington- 
Shepard instanton {N ^ oo, pi — p fixed), multiinstanton {N finite, pi finite) 
and dyon (pj = p ^ oo, N —>■ oo). 

When the radius of the loop R is much larger than the core of the solution 
(i.e. 3> p for (multi)instantons or b/27i for dyon), the Wilson loop is 

W{Cr) = exp(^iT3Ti^^^ = exp{iT3 ■ flux), (112) 

where Wr = f^||a?|=i,'. \x\ = r. 

Now for (multi) instanton one has for i? p 



Tir \r=R — ^ 

W 1 + E 



(113) 



B? + {xA-rib)'^ 



In nonsingular gauge one would obtain for (multi) instanton the flux 27r [68], 
in all gauges one has, that 

W{Cr) — 1, {multi)instantons. (114) 

Consider now the case of dyons, which amounts to tending p„ — > oo in Wr 
in (112). 

One can use the form (90) to obtain for dyon 

^ = -1; //«x = -7r, W{Cr) = -1. (115) 

It is amusing to consider also the intermediate case of so-called r-mono- 
poles [69], when p„ oo, but N is fixed, so the length of the chain L — Nb 
is finite. One can use (88) to find two limiting cases: 

i?>L, ^ = -2, flux ^-27r; W{Cr) = 1; (116) 
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R<^L, -^ = -1, flux^-n; W{Cr) ^ -1. (117) 

Thus only r-monopoles long enough, i.e. almost dyons, may ensure non- 
trivial Wilson loop, W{Cr) ^ 1. 

To connect flux values (114)-(117) to conflnement one can use model 
consideration of stochastic distribution of fluxes in the dilute gas, as was done 
in [68,69]. More generally the picture of stochastic fluxes was formulated in 
the model of stochastic confinement [70] and checked on the lattice in [71] . 

We shall come back to the model of stochastic conflnement in the next 
Section, and now shall use simple arguments from [68,69]. 

Indeed, consider a thin 3d layer above and below the Wilson loop, of 
thickness / -C R, and assume that it is flUed with gas of (multi)instantons 
or dyons. If 3d density of the gas is u, so that average number of objects is 
n — uSl, then Poisson distribution gives probability of having n objects in 
the layer around the plane of Wilson loop is 

w(n) = e"*^^. (118) 
nl 

If contribution to the Wilson loop of one object is A, A = +1 and —1 for 
instantons and dyons respectively, then the total contribution is 

< W{Cr) >= Yl e-"^^ = e-"(^-^) = e-"^, (119) 

where 

a = (l-A)z//. (120) 

Thus for instantons (A = 1) one obtains zero string tension in agreement 
with other calculations [55], while for dyons (A = — 1), conflnement is present 
according to the model. Let us stress, which features are important for this 
conclusion of conflnement: 

i) the flux =7r, so that W = —1 for one dyon, 

ii) stochastic distribution of fluxes, which enables us to use Poisson (or 
similar) distribution, 

iii) existence of flnite thickness, i.e. of flnite screening length, so that objects 
more distant than I completely screen each other and do not contribute to 
the Wilson loop. 
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Notice that point iii) is necessary for area law, otherwise (for / large, e.g. 
/ = R) one obtains a growing with R, i.e. super confinement. 

The same reasoning is applicable to torons and twisted instantons [66]. 
Indeed from (111) one can see, that their elementary flux is equal to (— tt). 
Therefore if one divides all volume into a set of twisted cubic cells, and 
ensures stochasticity of fluxes in the cells one should have the same result 
(119) with confinement present. 

Torons [72] and twisted instantons [66] have been studied from the point 
of view of confinement both analytically [72] and on the lattice [66]. For 
torons the requirement of stochasticity is difficult to implement, since bound- 
ary conditions of adjacent cells should ensure continuity of A^{x), and this 
introduces ordering in the fluxes, and confinement may be lost. In case of 
twisted instantons with Q = 1/2 [66] the field is essentially nonzero around 
the center of instanton, and b.c. are much less essential. The authors of [66] 
note a possibility of nonzero extrapolated value for the string tension when 
the box size is increased beyond 1.2. fm. It is still unclear what would be 
result when the twisted b.c. are imposed only on the internal boundaries. 

6.4. Below we concentrate on dyons, as most probable candidates for 

classical confining configurations. One must study properties of dyon gas 
and show that interaction in this gas is weak enough to ensure validity of 
the dilute gas approximation. As it is usual, one assumes the superposition 
ansatz 

N+ N- 

^. = E^^^^H^) + EV'^(^)' (121) 

i=l i=l 

where N^,N_ are numbers of dyons and antidyons respectively. To make 
the QCD vacuum 0(4) invariant, one should take any direction of the dyonic 
line, characterized by the unit vector uj^^ and position vector R^^\ so that 

Af(x)^nt(LA),{r,t)n,. (122) 

Here Qi is the color orientation matrix and L is the 0(4) (Lorentz) rotation 
matrix, while r and t are 

r^[{x- i?«)2 - {{x - i?«)^a;»)T/', (123) 

t= (a;-i?«)^a;«. (124) 
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It is now nontrivial, in which gauge take solution in (122), e.g. one may 
take singular gauge solution (91) or time-independent one (95), (96). The 
sum (121) is not obtained by gauge transformation from one case to another. 
Indeed it appears, that the form (95), (96) is not suitable, since the action for 
the sum (121) in this case diverges (see [67] for details). 

The form (91) is adoptable in this sense and wc consider it in more detail. 

Since solutions fall off fast enough, cf. Eq. (93), the interaction between 
dyons, defined as Sint, 

SiA)^ Y: S,iA^^) + Sir,t (125) 

i=l 

is Coulomb-like at large distances, e.g. for parallel dyon lines one has for 
two dyons 

where T — Nb is the length of dyon lines. A similar estimate can be obtained 
for nonpar allel lines. 

As was discussed in the previous section, the crucial point for the appear- 
ing of the area law of Wilson loop is the phenomenon of screening. To check 
this, let us consider the field of tightly correlated pair dd. When distance 
between d and d is zero, the resulting vector potential is obtained using the 
superposition ansatz (121) and d vector potentials (91)-(92 ) for d and the 
corresponding one for d, which differs from (91) by the sign of the last term 
and the total sign in (92): 

- 2 1 shr 

Aaidd) = -CaikUki- - cthr + — -), (127) 

g r chr — cost 

A4a{dd) = 0. (128) 

At long range one has 

Ai,{dd)^2eaik- + 0{e-'). (129) 
gr 

Calculation of B^a amounts to insertion in (98) f — which immediately 
yields: 

Bia{dd) = 0{e-'-), Eia = 0{e-'-). (130) 
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Thus fields of d and d completely screen each other at large distances; note, 
that this is purely nonabelian effect, since the cancellation is due to the 
quadratic term in / in (98). 

Now take distance between d and d equal to p — R'^^^ — R^'^^ and distance 
between observation point x and center of d J equal to r, f=x— M!^1±E^-^ 
assume that r ^ p. Then the field of dd (averaged over direction of p) is of 
the order 

Bk,Ek = 0{p^/r^). (131) 

Hence contribution of the distant correlated pair of dd is unessential, and 
indeed in calculation of the Wilson loop one can take into account the dis- 
tances to the plane of the loop smaller than the correlation length Z, which 
is actually the screening length. 

From the dimensional arguments - we have the only parameter in our 
Coulomb-like system - the average distance between the nearest neighbors 
therefore one has 

I = cl/-^/^ (132) 

where c is some numerical constant, and u is the 3d density of the dyon gas. 
Hence one expects the string tension in the dyonic gas to be of the order 

a = a?l'^. (133) 

Numerical calculations of < W{C^ > for the dyonic gas have been done in 
[73], but the density used was still much below that which is necessary for 
the observation of the screening; work is now in progress. 

Summarizing this Section, let us discuss perspectives of the classical so- 
lutions reported above as candidates for confining configurations. Only two 
solutions, dyons and twisted instantons, yield the suitable flux, equal to tt, 
through the Wilson loop, therefore we discuss them separately. Dyons can be 
represented as a coherent chain of instantons of large radius and correlated 
orientation of color fleld. When one goes from the instanton gas to these 
coherent chains , the action changes a little, but the enthropy decreases sig- 
nificantly, and possbly confinement occurs. To estimate the advantage or dis- 
advantage of dyonic configurations from the point of view of the minimum of 
the vacuum free energy, one must perform complicated computations, which 
are planned in the nearest future. 
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As to the twisted instantons, they require an internal lattice structure in 
the vacuum, which may violate the Lorentz invariance in some field correla- 
tors. 

After all, the problem is solved, as in dyonic case, by the calculation of the 
free energy of the vacuum: in the nature there should come into existence that 
vacuum structure, which ensures the minimal free energy. Lattice Monte- 
Carlo calculations satisfy this principle of minimal free energy (up to the finite 
size effects) and predict the confining vacuum with special nonperturbative 
configurations responsible for confinement. 

It is possible that those configurations are dyons, it is probable that they 
are not at all classical, but it is not excluded, that there exist unknown 
classical solutions, which ensure confinement after all. 



7 Topology and stochasticity of confinement 

In the previous chapter we have used the stochasticity of fluxes to obtain 

area law for dyons (magnetic monopoles). We start this chapter giving more 
rigorous treatment of this stochasticity and comparing it to lattice data . 

For abelian theory magnetic flux through the loop C is defined unam- 
bigously through the Wilson loop 

W{C) = exp ie j> A^dXfj, = exp ie / Hda, (134) 
J C J s 

and the magnetic flux is 

^i^e I Hda. (135) 

J s 

For SU(A'") theory the flux can be deflned analogously [70] (we omit the word 
"magnetic", since it depends on the orientation of the loop). 

Consider eigenvalues of the Wilson operator (note absence of trace in its 
definition) : 

U{C) = Pexp ig j> A/^dx^ = expiQ!(C). (136) 

The eigenvalues of the unitary operator U{C) are equal to expiQ;(C), where 
a{C) is a diagonal matrix N x N, depending on A^. 

In electrodynamics an{Ci2) are additive for the contour C12, consisting 
of two closed contours Ci and C2: 

«n(Ci2) =an(Ci)+a„(C2). (137) 
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In SU(A^) theory this is not so in general Consider now the spectral density 
Pc(Q!), i-e. averaged with the weight exp(— 5'o(A)) the probabihty of the flux 
a{C): 

^c(«) = JDA,e-si^) ' ^^^^^ 

where So{A) is the standard action of the SU{Nc) theory. 

Now any averaged Wilson operator over contour C, and also those for 
contour C", i.e. contour C, followed n times, can be calculated with the help 
oi pc{a): 

< W{C'') >= r dae"'''p,{a). (139) 



Assume, that there is confinement in the system, i.e. the area law holds both 
for contour C and C": 

< iy(C") >= exp(-A;„5). (140) 
Then the following equahty holds [70,71]: 

/TT PIT 
dai... / danPsl{ai)-Ps2io'n)S277{oi - ai - ... - q;„), (141) 
-TT J —IT 

where the contour C with the area S is made of contours Cj splitting the 
area S into pieces Si. The proof [70,71] can be done in both directions: from 
(139), (140) to (141) and back. Sometimes this statement is formulated as a 
theorem [70]: The necessary and sufficient condition 
of confinement is the additivity of random fluxes. 

The randomness is seen in (141), which has the form of convolution as it 
should be for the product of probabilities for independent events.The addi- 
tivity is evident from the argument of the 5-function in (141). 

The density Pc(ct) was measured in lattice calculations [71], and it was 
found, that Pc{<y) indeed satisfies (141) and approximately coincides with the 
pf^^(a) - density for d = 2 chromodynamics, if one renormalizes properly 
the charge. For the d — 2 case Pc{a) is also known explicitly and satisfies 
(141) exactly [71]. In that case confinement exists for trivial reasons. 

Let us now consider the nonabelian Stokes theorem [11], which for the 
operator (136) looks like 

U{C) = Pexpig j da^y{u)F^^{u,Xo), (142) 
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and take into account, that under gauge transformation V{x) it transforms 
as 



U{C) ^ V+{xo)U{C)V{xo), F^,{u,xo) ^ V+{xo)F^,V{xo). (143) 



Since U{C) can be brought to the form U{C) = expia with a diagonal by 
some unitary transformation, one can deduce, that this is some gauge trans- 
formation V{xq), and, moreover, this is the same, which makes / d(7^,^F^,^(u, xq) 
diagonaL Thus one can define the flux /i similarly to (135): 



Note, that dependence on xq presents in U{C), cancels in W{C) — trU{C). 
The additivity of fluxes is seen in (144) explicitly. 

Now consider the statistical independence of fluxes, which obtain, when 
one divides the surface S into pieces Si,...,Sn- Using the cluster expansion 
theorem [12] and discussion in Section 3, one can conclude, that necessary 
and sufficient condition for this is the finite correlation length Tg, which ap- 
pears in correlation functions (cumulants) <^ F{l)...F{n) In this case, 
when pieces Sk, k = 1, ...,n, are all much larger in size than Tg, then different 
pieces become statistical independent. Thus our consideration in the frame- 
work of the field correlators in Section 3 is in clear agreement with the idea 
of stochastic confinement [70,71]. The MVC in addition contains the quanti- 
tative method to calculate all observables in terms of given local correlators, 
which is absent in the stochastic confinement model [70,71]. 

The appearance of new physical quantity - Tg and cumulants is a fur- 
ther development of the idea of stochastic vacuum, which gives an exact 
quantitative characteristics of randomness. 

When size of contours is of the order of Tg, the fluxes are no more random, 
and area law at such distances disappears - there is no area law at small 
distances, as is explained in Section 3. 

The lattice measurements in [71] show that Pc{o:) is strongly peaked 
around a = for small contours, when there is no area law, which corre- 
sponds to the perturbative regime. 

Instead for large contours the measured Pc{<y) are rather isotropic. 

This fact one can compare with our result for fluxes for instantons and 
dyons. From our deflnition of fluxes (136) dyon in the plane of contour C 




(144) 
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corresponds to (cf.Eq.(112)) 




(145) 



and dyon with the center off the plane has smaller eigenvalues cc^. It is clear, 
that instantons with flux zero cannot bring about isotropic distribution of 
fluxes, while having maximal flux (±7r) are most effective in creating the 
isotropic pc{ot)-i when one integrates over all dyons in the layer above and 
below the plane. 

It is instructive now to study the question of fluxes for adjoint Wilson 
loop and in general for Wilson loops of higher representations. 

One can keep the dcflnition (136) also in this case, but and a.{C) 
should be expressed through generators of given representation: 



(146) 



Thus o{C) for the Wilson loop in adjoint representation is a matrix 
{N^ - 1) X (iV2 - 1). E.g. for SU(2) (T")f,c = |e„6,. To understand how 
stochastic wacuum model works for the adjoint representation, let us take 
as an example the flux of one dyon and calculate oiadjiC). Repeating our 
discussion, preceding eq.(112) for large loops in the (12) plane, one concludes, 
that again only color index a — 3 contribute, and one has 



aadj{C) = TT • diag{T^)lim(^^^ 



Since the last factor for dyon is lim^^ — —1, and diag{T^) 
one has flnally 



/ 1 

V 



— TT 



TT 



(147) 



(148) 







Thus our conclusion of the elementary flux equal to n holds true also 
in the adjoint representation (and all higher representations), which gives 
an argument for confinement of adjoint charges on the same ground as for 
fundamental charges. 
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One has 

< Wadj >= j^Tzri^^-dje'^"'^''^ = (149) 

as well as < Wfund >— — 1- 

So far we have discussed stochasticity of the vacuum from the point of 
view of fluxes and conclude, that it shows up as random distribution of 
fluxes. In Section 3 the vacuum stochasticity was formulated in the language 
of field correlators. Through the AP method one can connect the latter with 
the distribution of AP magnetic monopole currents, (in the U{1) theory an 
exact connection holds even without AP). One may wonder, why magnetic 
monopoles or dyons are needed to maintain the stochastic picture of the 
vacuum? 

To answer this question we start with the Abelian theory. Without mag- 
netic monopoles Bianchi identities divH — are operating, requiring, that 
all magnetic field strength lines are closed. 

This introduces strong ordering in the distribution of magnetic field, and 
no stochastic picture emerges. 

As a result, confinement is not present in the system, as can be seen from 
eq.(22). In presence of magnetic monopoles the magnetic hues can start 
and end at any place, where monopole is present, and one can have a really 
stochastic distribution. As we discussed it in Section 3, the nonabelian dyna- 
mics can mimick the effect of monopoles due to triple correlators < EiEjBk > 
and ensure in this way the stochastic distribution of fields. 

So magnetic monopoles in Abelian theory, dyons in gluodynamics create 
disorder in the system. 

The same situation occurs in other spin and lattice systems, e.g. in the 
planar Heisenberg model the Berezinsky-Costerlitz-Thouless vortices cre- 
ate disorder and and master the phase transition into the high-temperature 
phase [74] (for details see also [1]). 

All that is a manifestation of a general principle [2] : 

Topologically nontrivial field configurations are responsible for creation 
of disorder and they drive the phase transition order-disorder. 

From the QCD point of view the "ordered phase" is the perturbative vacu- 
um of QCD with long distance correlations {Di (x) ~ ^) and flux distribution 
Pc{c() centered at zero, while the "disordered phase" is the real QCD vacuum 
with short correlation length Tg and with random fluxes; there is no real 
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phase transition in continuum: the phases coexist on two different scales 
of distances (or moments). In the lattice version of U{1) theory there are 
indeed two phases: weak coupling phase, corresponding to the usual QED, 
and strong coupling phase with magnetic monopoles - lattice artefacts - 
driving the phase transition. 

The dyon is a continuum example of topologically nontrivial configu- 
ration. In singular ('tHooft's) gauge dyon has multiinstanton topological 
number, proportional to its length. 

Dyon saturates the triple correlator < EiEjBh > and may be a source of 
randomness of field distribution. 

However the ultimate answer to the question about the nature of confin- 
ing configurations is still missing. The analysis of the dyonic vacuum as a 
model of the QCD vacuum is not yet completed, and it is possible that the 
topologically nontrivial configurations responsible for confinement are dyons, 
or some other not known solutions or else purely quantum fluctuations. 

We conclude this Section with the discussion of a possible connection 
between confinement and the Anderson localization [75]. 

At the base of the similarity between these two phenomena lies the field 
stochasticity in the vacuum (medium), where quark (electron) propagates. 
The similarity however ends up just here. 

Namely, for an electron one can discuss its individual Green's function, 
which always (for any density of defects) decays exponentially with distance. 
[There exists, however, a special correlator, e.g. the direct current conduc- 
tivity adc, which vanishes for localized states ( for large density of localized 
defects [76]) and is nonzero for the delocalized states]. 

In the case of a quark in the confining vacuum its Green's function (more 
precise: the gauge-invariant Green's function of the qq system averaged over 
vacuum configurations) corresponds to the linear potential, i.e. it behaves 
as G{r) ~ exp(— r^/^), where r is the distance between q and q. Thus the 
quark Green's function decays faster than any exponent, in contrast to the 
Green's function of an electron in the medium, always decaying exponentially. 
This property of vacuum Green's functions was coined by the author [77] the 
super localization. If the average potential V, acting on a quark, had been 
finite, then quark at some high energy could be freed and get to the detector. 

The essence of the superlocalization is exactly the fact, that the averaged 
potential V grows with distance without limits and therefore quarks are 
confined at any energy -this is the absolute confinement . 
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It is interesting to follow the mechanism, how the unbounded growth of 
V occurs. To this end consider, as we did at the beginning of Section 3, a 
nonrelativistic quark, moving in the x, t plane, while the heavy antiquark is 
fixed at the origin. 

According to the quantum mechanical textbooks [78] the quark Green's 
function is proportional to the phase integral, and using the Fock-Schwinger 
gauge, one can write 



1-^ 



G{X, T) =< exp ig J ^4(2;, t)dt >^ 

T T J'^ 

[ dt [ dt' [ du I du' < Ei{u,t)Ei{u\t') >K (150) 

Jo JQ Jo Jo 



1 - VT 



Stochasticity of vacuum fields implies the finite correlation length Tg for the 
correlator < Ei{u,t)Ei{u' ,t') >, i.e. according to (16) one has 

< Ei{u, t)Ei{u', t') >= D{u -u',t- t') + ... (151) 

and for large T and X we obtain 

VKC(mst\X\, \X\^oo (152) 

Thus the linear growth of is a consequence of random distribution of 
field strength E{u, t) and of the fact, that y is a result of the averaging of 
vector-i)otoiitials A^, which are connected to F^^, by an additional integral. 
This extra integration causes the linear growth of V , and from the physical 
point of view this means the accumulation of fluctuations of the field F^,, on 
the whole distance X from the quark to the antiquark. 

This is the essence of the superlocalization phenomenon, which still has 
no analogue in the physics of condensed matter. 



8 Conclusions 

In this rewiew we have looked at conflnement from different sides and de- 
scribed the mechanism of this phenomenon in the language of field corre- 
lators, using more phenomenological language of dual superconductivity - 
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effective classical equations of the Ginzburg -Landau type, in the language 
of the stochastic flux distributions and finally we have studied classical con- 
figurations which may be responsible for confinement. 

All the way long we have stressed that confinement is the string forma- 
tion between color charges, the string mostly consisting of the longitudial 
colorelectric field. 

Let us try now to combine different descriptions of confinement, given 
above in the review, and show in a simple example how the string looks like. 

To this end we use the simple picture of the nonrelativistic quark and 
the heavy antiquark at the distance X between them, discussed at the end 
of the last Section. From the point of view of field correlators, confinement 
- the string formation - is the consequence of the fact, that there exists the 
correlation length Tg, such that fields inside this length are coherent and those 
outside of this length are random. This is shown in Fig. 10.1 (left part), where 
the strip of size Tg is indicated in the plane (14) (one could take instead any 
other plane, i.e. the plane (12) or (13)). Inside this strip the field is directed 
mostly in the same way (i.e. as on the string axis), but outside of the strip 
directions are random. The correlation length Tg characterizes the string 
thickness (if the plane (12) or (14) is chosen). Thus Tg plays the double role, 
it gives the coherence lengh, where the string is created, and beyond which - 
stochastic vacuum field, existing also before quark and antiquark have been 
inserted in the vacuum. 

Let us look at the same construction from the point of view of the dual 
superconductivity. Then one obtains the picture shown in Fig. 10. 2. Here an 
arrow depicts the monopolc current j^, which is caused by the colorelectric 
field Ex of the string in accordance with the dual London's equation rotj = 
m^E. The effect of this current is the squeezing of the string field, which 
prevents the field flux lines from diverging into the space, and as a result E^ 
decays exponentially away from the string axis Ox like exp{—m\/y^ -\- z'^). 
Thus m defines the string thickness and one can conclude that m ~ ^/Tg. 
And indeed Eq.(22) confirms this conclusion. 

Let us turn now to the fiux distribution and to the stochastic vacuum 
model. Fig. 10.3. In this case the strip, corresponding to the string in Fig. 
10.1 can be divided into pieces Si, S2, ... of the size (P, such that fluxes inside 
each piece are coherent and equal to e.g. ±7r for the case of dyons, but two 
neighboring pieces are noncoherent - their fluxes are mutually random. The 
string thickness is now built up due to the size d of the piece Sm containing 
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a coherent flux. If tlie surface Sn is penetrated by a monopole or a dyon, 
tlien d coincides with the monopole or dyon size. 

To clarify this point let us find the minimal size of the loop R where the 
dyon flux is equal to the asymptotic value of (— tt). To this end we use (115) 
and insert there (90), and obtain the result, that for R ^ 6/27r = 7"^ the 
flux is equal to (— tt) with exponential accuracy; hence the size of the dyon 
flux is equal to 7"^ and this should be the size d of the piece Sn- 

Prom the point of view of field correlators d should coincide with Tg, 
therefore the string thickness is of the order of the typical dyon (or monopole) 
size (or an other classical solutions). 

Hence, all our pictures represented in Fig. 10.1- 10.3 can be combined 
in one generalized mechanism of the string formation, which is based on the 
existence of coherent field domains of the size Tg, and beyond that size the 
fields are independent and random. 

A question arises: who manages this structure of QCD vacuum and why 
in the case of QCD and gluodynamics the vacuum is made up this way, while 
in the case of QED and the Weinberg-Salam theory , the nonperturbative 
configurations are probably suppressed and the vacuum structure is differ- 
ent. To be able to answer this question is also to answer the question of 
phase transition mechanism for the tempcratiu^c dcconfinemcnt, which was 
observed on the lattice [79]. This topic requires a separate review paper, 
since the amount of information accumulated here by now is very large. We 
shall confine ourselves to only few remaks on this point. 

Firstly, the density of the (nonperturbative) vacuum energy one can con- 
nect with the help of the scale anomaly theorem with the magnitude of the 
nonperturbative gluonic condensate [8] 

e = +^ < i^;.(0)F;.(0) > (153) 

For small ctg the function (3{as) is negative - in contrast to QED, and if 
it keeps the sign in all effective region of as, then one can deduce that the 
nonperturbative vacuum shift (153) is advantageous, since it diminishes the 
vacuum energy ( and also the free energy for small temperatures). 

This conclusion can be considered as an intuitive idea why the nonper- 
turbative vacuum in QCD is advantageous and comes into existance, while 
QCD - not advantageous and is not realized. 
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Secondly, let us briefly discuss the phase transition with increasing tem- 
perature in QCD, referring the reader to lattice calculations [79] and original 
papers [81] for details, The main criterium which defines the vacuum struc- 
ture preferred at a given temperature, is the criterium of the minimum of the 
free energy (which is a corollary of the second law of thermodynamics). In 
the confining phase for T > the free energy consists of the term (153) and 
of the contribution of hadronic excitations (glucballs, mesons and baryons), 
which slowly grows up to T ~ 150MeV. Note that the gluonic condensate 
contains both colorelectric and colormagnetic fields, but only the first ones 
have to do with confinement in the proper sense of this term. 

The deconfinement phase, realized at T > T^, usually was identified as 
the phase with perturbative vacuum, where quarks and gluons in the lowest 
order in g are free [82]. However, from the point of view of the minimum 
of free energy it is advantageous to keep in the vacuum colormagnetic fields 
and the corresponding part of the condensate (153) since quarks in this case 
stay practically free, and one significantly gains in energy - around one half 
of the amount in (153). This is the "magnetic confinement" phase [81]. 
Calculations in [81] yield Tc in good agreement with lattice data for number 
of flavours Uf — 0, 2, 4. 

The main prediction of the "magnetic conflnement" is the area law for 
T > Tc of the spacial Wilson loops [83], and the phenomenon of the "hadronic 
screening lengths" , i.e. the existence of hadronic spectra for T > Tc in the 
Green's functions with evolution along space directions [81], which agrees 
well with the lattice measurements [79] . 

Hence the picture of the phase transition [81] into the "magnetic conflne- 
ment" phase, supported by computations, seems to be well founded. In this 
picture at T > Tc disappear colorelectric correlators , more explicitly, the 
correlators of the type of of D{x), contributing to the string tension. 

What happens then with effective or real magnetic monopoles and dyons? 
In the AP method at T > Tc the monopole density strongly decreases [35], 
which can be understood as an active annihilation or a close pairing of 
monopoles and antimonopoles. The same can be said about pairs of dyons 
and antidyons. Thus the deconflnement phase of color charges can be asso- 
ciated with the conflnement phase of monopoles (or dyons). However, the 
"magnetic conflnement" phenomenon imposes deflnite requirements on the 
vacuum structure at T > Tc . E.g. there should exist magnetic monopole 
(dyonic) currents along the 4-th (euclidean time) axis -i.e. static or periodic 
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monopoles (dyons). Such currents can confine in spacial planes (what is ob- 
served on lattices [83]), but do not participate in the usual confinement (i.e. 
in the temporal planes (i4),i = 1,2,3). These points will be elucidated in 
separate publications. 

We also had no space to discuss an important question of the connec- 
tion between confinement and spontaneous breaking of chiral symmetry and 
Ua{^), where magnetic monopoles (dyons) can play an important role [84]. 

We always held above that confinement is the property not only of QCD 
( with quarks present in the vacuum), but also of gluodynamics (without 
quarks). This conclusion follows from numerous lattice data (see e.g. [10]), 
and also from computations supporting the dual Meissner effect as a basis of 
confinement discussed in the review, where quarks do not play an important 
role. 

On the other hand there arc not any calculation or experimental data 
which support a key role of quarks in the confinement mechanism. For this 
reason we have not discussed above the model of V.N.Gribov (interesting by 
itself) and the reader is referred to original papers [85]. 
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and to K.A.Ter-Martirosyan for constant support and discussions. The fi- 
nancial support of RFFR, grant 95-02-05436, and of INTAS, grant 93-79 is 
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Figure captions 



Fig. 1. Potential between static quarks in the triplet representation of 
SU(3), computed in [10] on the lattice 32^. The solid line - the fit of the 
form ^ + aR + const. The potential and distance R are measured in lattice 
units a, equal to 0.055 fm for {3 = 6.9. Dynamical quarks are absent. 

Fig. 2. The same potential as in Fig. 1, but with dynamical quarks taken 
into account in two versions (staggered fermions - upper part, a ^ 0, 11 fm, 
Wilson fermions - lower part, a ~ 0, 16 fm); calculations of Heller et al (sec- 
ond entry of [10]). 

Fig. 3. The same potential as in Fig. 1, but for quarks in the sextet (a) 
and octet (b) representation. Broken line - the triplet potential of Fig. 1 
multiplied by the ratio of Casimir operators, equal to 2.25 for the octet and 
2.5 for the sextet. 

Fig.4. Correlators Du{x) — D + Di + x^^- (lower set of points) and 
D±{x) — D + Di (upper set of points) as functions of distance x. Crosses 
correspond to d = 5.8, diamonds - to = 5.9 and squares - to f3 = 6.0. 
Solid lines are best fits in the form of independent exponents for D{x) and 
Di{x). Computations from [13]. 

Fig. 5. Distribution of the parallel colorelectric field Eu as a function 
of distance x± to the string axis. Measurements from [19] are made at dif- 
ferent distances Xn from the string end: filled circle - Xu = 3a, diamond 
-Xii — 5a and triangle -Xu = 7a (a is the lattice unit ) . Solid line - Gaus- 
sian fit, dashed line - calculation in [19] with the help of D{x) taken from [13]. 

Fig. 6. The same distribution En, as in Fig. 5, for SU(2) gluodynamics 
measured in [24] on the lattice 24"^ for (3 = 2.7, Xn = 5a, (the length of 
the string is 10a) - empty circles, compared to the Bu distribution in the 
Abrikosov string - solid line. The growth od the solid fine at small Xt, is 
unphysical and is due to violation of approximations made in (47). 

Fig. 7. Lattice measurements [25] of the penetration length 6 = X and 
the coherence length ^ as functions f5 = 2Nc/g'^ for AP configurations in 
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SU(2) gluodynamics (upper two figures) and for SU(3) (tlie lower figure). 
Tlie values of 6 and ^ are defined by comparison of field distribution and AP 
monopole currents with the solution of the Ginzburg-Landau equations. 

Fig. 8. The string tension measured in [43] for all AP configurations 
(squares) and separately for AP monopoles (empty circles) and "photons" 
(filled circles) as functions of P = A/g^ in SU(2) gluodynamics. 

Fig. 9. Effective potential of the AP monopole field (p, defined according 
to [48], measured in [26] for two values of (3 in SU(2) gluodynamics, corre- 
sponding to confinement (Fig.9a) and deconfinement (Fig.Qb). 

Fig. 10. The picture of string formation between a nonrelativistic quark 
and a heavy antiquark, is illustrated in three different approaches, discussed 
in the text. 

(1) - in the formalism of field correlators 

(2) - in the formalism of dual superconductivity 

(3) - in the picture of the stochastic flux distribution. 



57 



References 

[1] Simonov Yu.A. Models of confinement, Proceedings of the XXII LINP 
Winter School (L:LINP,1985) 

[2] Polyakov A.M. Nucl. Phys. B120, 429, 1977 
Polyakov A.M. Phys. Lett. 59, 80,82, 1975 

[3] 't Hooft G. in High Energy Physics (Ed. Zichichi A.) (Bologna: Editrice 
Compositori, 1976) 

Mandelstam S. Phys. Lett. B53, 476, 1975 

[4] 't Hooft G. Nucl. Phys. B190,455, 1981 

[5] Dosch H.G. Phys.Lett. B190, 177, 1987 

Dosch H.G., Simonov Yu.A. Phys.Lett. B205, 339, 1988 
Simonov Yu.A. Nucl. Phys. B307, 512, 1988 
Simonov Yu.A. Yad.Fiz. 54, 192, 1991 

[6] Wilson K.G. Phys.Rev. DIO, 2445, 1974 

[7] Greensite J., Halpern M.B. Phys.Rev. D27, 2545, 1983 

[8] 't Hooft G. Nucl. Phys. B72, 461, 1974 

[9] Ellis R.K. QCD at TAS'94, Fermilab - Conf-94/410-T (Fermilab, 94) 
AltareUi G. in: QCD - 20 years later - Eds. Kastrup H.A., Zerwas P. 
(World Scientific, 1993) 

[10] Ford l.J. et al Phys. Lett. B208, 286, 1988 Heller U.M. et al. Phys. Lett. 
B 335, 71, 1994 

[11] Arefieva I.Ya. TMP 43, 111, 1980 

Simonov Yu.A. Yad.Fiz. 50, 213, 1989 

[12] Van Kampen N.G. Phys.Rep. 24C, 171, 1976 
Physica 74, 215, 1974 
Simonov Yu.A. Yad.Fiz. 48 , 1381, 1988, 
Yad.Fiz. 50 , 213, 1989 

[13] Di Giacomo A., Panagopoulos H. Phys.Lett. B285, 133 , 1992 



58 



[14] De Grand T.A., Toussaint D. Phys. Rev. D22, 2478, 1980 
D24, 466, 1981 

Bhanot G. Phys. Rev. D24, 461, 1981 



[15 
[16 

[17; 

[18 
[19 

[20 

[21 
[22 

[23 

[24 
[25 

[26 
[27 
[28 
[29 
[30 
[31 



Simonov Yu.A. Yad.Fiz. 50 , 213, 1989, 
Durand L., Mendel E. PPhys. Lett. 85B, 241, 1979 
Simonov Yu.A. Yad.Fiz. 50, 500, 1989 
Simonov Yu.A. Yad.Fiz. 57, 1491, 1994 

Del Debbio L., Di Giacomo A., Simonov Yu.A. Phys.Lett. B 332, 111 , 
1994 

Lee T.D. Particle Physics and Introduction to Field Theory.- 1981. Har- 
wood. New York. Chapter 17. 

Alder S., Piran T. Rev.Mod. Phys. 56, 1, 1984 

Nielsen H.B., Olesen P. Nucl. Phys. B57, 367, 1973 
B61, 45, 1973 

Lifshitz E.M.,Pitaevski L.P. Statical Physics, 4.2. Chapter 5, Moscow, 
Nauka, 1973 

Cea P., Cosmai L., Nucl. Phys. B (Proc. Suppl.) 42, 225, 1995 

Matsubara Y.,Ejiri.S., Suzuki T. Nucl. Phys. (Proc. Suppl.) B34, 176 , 
1994 



Chernodub M.N., Polikarpov M.I., Veselov A.L |hep-lat/9512"030| 



Polikarpov M.I., Uspekhi Fiz.Nauk, 165, 627, 1995 

Tinkham M. Introduction to Superconductivity, (Kreiger, USA, 1980) 

ref. [23], §26 

Weinberg S. Progr. of Theor. Phys. Suppl. 86, 43 , 1986 
Molodtsov S.V., Simonov Yu.A. Pis'ma ZhETF 60, 230, 1994 

59 



[32 
[33 
[34 
[35 
[36 

[37; 

[38 

[39 
[40 
[41 

[42 
[43 

[44; 
[45 
[46 



Simonov Yu.A. Yad.Fiz. 58, 113, 1995 

Baker M. Ball J.S., Zachariasen F. Phys. Rev. D37, 1036, 1988 
Baker M. et al. Phys. Rev. D51, 1968, 1995 

Kronfeld A.S. Schierholz G., Wiese U. Nucl. Phys. B293, 461, 1987 

Chernodub M.N., Polikarpov M.I., Veselov A.I. Phys. Lett. B342, 303, 
1995 

Suzuki T. Nucl. Phys. (Proc.Suppl.) B30, 176 , 1993 

Suzuki T. et al. Phys. Lett. B347, 347 , 1995 

Stack D. et al. D50, 3395, 1994 

Ejiri S. et all. Nucl. Phys. B42, 481, 1995 

Suzuki T. in Lattice - 95 

Miyamura O., Origuchi S. |hep-lat/9508015 



Ken Yee. Phys. Lett. B347 367 1995 

Suganuma H. et al [hep-ph 9506366| , in Lattice - 95 

Brown D.A., Haymaker R.W., Singh V. Phys. Lett. B306, 115, 1993 

Matsubara Y. et al. Nucl. Phys. (Proc.Suppl.) B34, 176 , 1994 
Shiba H., Suzuki T. Nucl. Phys. (Proc.Suppl.) B34, 182 , 1994 



Haymaker R.W., Lectures at E.Fermi School, Varenna 1995 
Cea P., Cosmai L. in Lattice - 95 

Schichter C, Bah G.S. , Schilling K. Nucl. Phys. B (Proc. Suppl.) 42, 
273, 1995 



[47] Ivanenko T.L., Polikarpov M.I., Pochinski A.V. Phys. Lett B252, 631, 
1990 

ibid B302, 458, 1993 

Di Giacomo A. plenary talk at Lattice - 95 



60 



[48] Froehlich J., Marchetti P.A. Comm. Math. Phys. 112, 343, 1987 

[49] Polley L., Wiese U.J. Nucl. Phys. B356, 629, 1991 

Pohkarpov M.I., Polley L., Wiese U. Phys. Lett. B253, 212, 1991 

[50] Del Debbio L. Di Giacomo A., Paffuti G. Phys. Lett. B349, 513 , 1995 
Nucl. Phys. (Proc.Suppl.) B42, 231 , 1995 

[51] Del Debbio L. et al. Phys. Lett. B355, 255 , 1995 
Nucl. Phys. (Proc.Suppl.) B42, 234 , 1995 

[52] Chernodub M.N., Gubarev F.V. Pisma v ZhETF 62, 91, 1995 

[53] Hart A., Teper M. hep-lat 95 11016 

[54] Bornyakov V., Schierholz G. (in preparation) 

[55] Callen C.G., Dashen R., Gross D.J., Phys. Lett. B66, 375, 1977 
Veselov A.I., Polikarpov M.I., Pisma v ZhETF 45, 113, 1987 

[56] Campostrini M., Di Giacomo A., Phys. Lett. B225, 403, 1989 

[57] Jackiw R., Nohl C., Rebbi C. Phys. Rev. 15, 1642, 1977 

[58] Atiyah M.F., Kitchen N.J., Drinfeld V.G., Manin Yu.I. Phys. Lett. A65, 
185, 1978 

[59] Belavin A. A., Polyakov A.M., Schwartz A.S., Tyupkin Yu.S. Phys. Lett. 
B59, 85, 1975 

[60] Harrington B., Shepard H. Phys. Rev. D17, 2122, 1978 
[61] Rossi P. Phys. Rep. 86, 317, 1982 

[62] Prasad M.K., Sommerfield CM. Phys. Rev. Lett. 35,1760, 1975 

[63] Polyakov A.M., JEPT Lett., 20, 894, 1974 
't Hooft G. Nucl. Phys. B79,276, 1974 

[64] 't Hooft G. Nucl. Phys. B153,141, 1979 
Comm. Math. Phys. 81, 267, 1981 



61 



[65] Leutwyler H. Nucl. Phys. B179,129, 1981 

[66] Garcia-Perez M., Gonzalez-Arroyo A., Martinez P. Nucl. Phys. (Proc. 
Suppl.) B34,228, 1994 

[67] Simonov Yu.A. Dyons in QCD: confinement and chiral symmetry break- 
ing - Lectures at the Intern. E.Fermi School, Varenna, 1995 (in press) 

[68] Simonov Yu.A. Yad.Fiz. 50, 500, 1989 

[69] Simonov Yu.A. Yad.Fiz. 42, 557, 1985; 
preprint ITEP-156, (M.: ITEP, 1985 ) 

[70] Olesen P. Nucl. Phys. B200, 381, 1982; 

Ambjorn J. Olesen P., Peterson V. Nucl. Phys. B240, 553, 1984 
Ambjorn J. Olesen P. Nucl. Phys. B170, 60, 265, 1980 

[71] Makccnko Yu.M., Polikarpov M.I., Veselov A.I. Phys. Lett. B118, 133, 
1982 

Belova T.I., Makeenko Yu.M., Pohkarpov M.I., Veselov A.I. Nucl. Phys. 
B230, 473, 1984 

[72] Simonov Yu.A. Yad.Fiz. 46, 317, 1987 

[73] Veselov A.I., Martemianov B.V., Molodtsov S.V., Simonov Yu.A. Pis'ma 
V ZhETF , 62, 679, 1995 

[74] Berezinski V.L. ZhETF 59, 907, 1970; 
61 1144, 1971 

Kosterlitz J.M., Thouless D.J. J.Phys. C6,118, 1973 

[75] Anderson P.W. Phys. Rev. 109, 1492, 1958; 
Rev. Mod. Phys. 50, 191, 1978 

[76] Lee P. A., Ramakrishnan T.V. Rev. Mod. Phys. 57, 287, 1985 

Lifshitz I.M., Grcdcskul S.A., Pastur L.A. Introduction into theory of 
disordered systems (in russian) (Moscow, Nauka, 1982) 

[77] Simonov Yu.A. Preprint ITEP 68-87 (M.: ITEP, 1987) 



62 



[78] Feynman R., Hibbs A. Quantum Mechanics and Path Integrals Mc 
Graw-Hill, New York, 1965) 

[79] Karsch F. in "QCD: 20 years later" - Eds. Zerwas P.M. , Kastrup H.A., 
(World Scientific 1993) 

[80] Vainstein A.I., Zakharov V.I. Novikov V.A., Shifman M.A. Uspekhi Fiz. 
Nauk 136, 553, 1982 

[81] Simonov Yu.A. Pis'ma ZhETF 55, 605, 1992, Yad.Fiz., 58, 357, 1995 
Dosch H.G. , Pirner H.J., Simonov Yu.A. Phys.Lett. B349, 335, 1995 
Gubankova E.L., Simonov Yu.A. Phys. Lett. B360, 93, 1995 

[82] Kapusta J. Finite Temperarure Field Theory - Cambridge Univ. Press. 
1989 

[83] Manousakis E., Polonyi J. Phys. Rev. Lett. 58, 847, 1987 
Bah G.S. et al. Phys. Rev. Lett. 71, 3059, 1993 

[84] Gonzalez- Arroyo A., Simonov Yu. A. Nucl. Phys. B (in press) 
|hep-th 95060321 

[85] V.N.Gribov, Physica Scripta T15, 164, 1987; preprint LUTP 91-7 (Lund 
Univ., 1991) 



63 



